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Introduction 


In Unit 7 you studied a model of vibrating systems which is based upon the concept of 
a perfect spring and leads to the prediction of simple harmonic motion. This unit 
continues the theme of modelling vibration, introducing the three ideas of damping, 
forcing and resonance. | 


Damping refers to the inevitable tendency of vibrations to die away whenever a 
system is left to its own devices. This phenomenon was clearly apparent in the 
oscillations of the home-made apparatus of Unit 7 Section 1, but no attempt was made 
there to include damping in the model. In this unit we model damped vibrations by 
including a velocity-dependent resistive force in Newton’s second law. 


Forcing is what is needed to keep a system vibrating, in spite of its damping. This can 
arise from a periodic force. For example, a car engine vibrates due to the rotating 
crank even when the car is stationary. A mechanical system can also be forced to 
vibrate by the specific movement of some component of the system. Thus a car on a 
bumpy road is forced to vibrate, at a frequency not of its own choosing, by the bumps 
which its wheels encounter. Similarly, the tone-arm of a record player responds to the 
motion of the stylus in contact with a record. 


Resonance refers to the fact that the amplitude of forced vibrations can be much 
greater at some forcing frequencies than at others. These resonant oscillations can be 
destructive in mechanical systems. Often designers introduce springs and damping 
devices in order to eliminate resonant behaviour over the range of frequencies which a 
piece of machinery is likely to encounter. 


The mathematics required here is that which you studied in Unit 6. In particular, the 
mechanical systems to be considered are modelled by the linear differential equation 


mz+re+kx = kx, + Pcosuwt, 


where m, r, k, x-, P and w are constants. The case for which r = 0 and P= 0 was ~ 
discussed in Unit 7, and corresponds to simple harmonic motion. Damped vibrations 
are modelled by the addition of the term rz, where r is a positive constant. The 
inclusion of the term P coswt corresponds to forced vibrations. 


Study guide 


The first four sections of this unit should be studied in the order that they appear. The 
main ideas are in Sections 1 and 2. Section 1 uses Newton’s second law, with perfect 
spring forces and resistive forces which are linearly dependent on velocity, to model 
damped vibrations. A periodic applied force is added to this model in Section 2 in order 
to incorporate forced vibrations. This basic model is generalized in two respects in 
Section 3. First, damping forces are discussed which depend on the relative velocity 
between two components of a system. Secondly, we consider forced vibrations which 
arise from the periodic movement of a point of the system other than the particle 
whose motion is sought. 


The audio-tape is associated with Section 3. It will give you practice in using Newton’s 
second law to derive the equations of motion which model mechanical systems. The 
television programme is associated with Section 4 and provides a basis for evaluating 
the success of our models of damped and forced vibrations. It will be easier to follow if 
you have already studied the first three sections. 


Finally, Section 5 contains some further exercises on the whole unit. You may use these 
for additional practice, either during the current study week or when you are revising 
for the examination. 
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1 Damped vibrations 


1.1 Introduction 


In Unit 7 we modelled a variety of vibrating systems, and in each case the model led to 
the prediction of simple harmonic motion. But in the real world no system, left to its 
own devices, maintains the oscillations of constant amplitude that are predicted by 
such models; in reality the amplitude decreases with time. In this section we therefore 
introduce a new ingredient to model the observed behaviour of damped vibration. This 
ingredient is a velocity-dependent force. 


One particular system considered in Unit 7 was a home-made apparatus used in a 
simple experiment. The system consisted of a bag of coins suspended from an elastic 
string whose other end was attached to a fixed frame. The bag was pulled down below 
its position of rest and then released, after which it oscillated up and down between 
extreme positions on either side of its equilibrium position. The mathematical model of 
Unit 7, used in an attempt to describe this behaviour, predicts that when the 
displacement of the bag of coins is plotted against time, the resulting curve will be 
sinusoidal, as shown in Figure 1. 


Ree Ws eae ig ieee eae. 
A 
Se) Se See ee eS ee Ne | 
A 
| period 7 = 27/w 
O t 


Figure 1 Simple harmonic motion with amplitude A and period rT. 


Specifically, the model predicts that the equation of motion is % + w*x = wx, and that 
the displacement of the bag is given by x = x, + Acos(wt + ¢), where x, is the 
equilibrium displacement, A is the amplitude, w is the angular frequency and ¢ is the 
phase angle. 


height of bag 
above ground 


Bie ee ee are ee er 
height 


time 


Figure 2. The observed behaviour of the home-made apparatus of Unit 7 Section 1, 
showing damped vibration. 
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For this simple harmonic motion both the amplitude and the period t = 27/w are 
constant. How does this prediction tally with observed fact? The experimental results 
shown in Figure 3 of Unit 7 Section 1 supply the answer. The essential features of this 
observed motion are repeated in Figure 2 above. 


In one respect the prediction corresponds with what is observed: the motion does have 
a constant period. But from another point of view the prediction is clearly wrong: the 
motion is not one with a constant amplitude. In fact, the amplitude decreases steadily 
with time and eventually the motion ceases altogether. This observed decrease in 
amplitude is not represented in the simple harmonic model, and in the current section 
we intend to remedy this deficiency. 


1.2 Linear damping 


It is a fair guess that a major reason for the decreasing amplitude of the vibration 
represented in Figure 2 is the effect of air resistance. However, it is known that the 
repeated lengthening and shortening of a rubber band causes an energy loss due to 
internal friction, and this will also contribute to the decrease in amplitude. In the 
context of vibration, frictional effects are commonly referred to as damping and the 
corresponding motion is spoken of as a damped vibration. 


Damping can be due to a number of physical effects and there is a corresponding 
variety of mathematical models. In this unit we concentrate on the resistance 
experienced by a solid body moving through a liquid or gas. As you saw in Unit 4, this 
resistance is commonly modelled as a force opposed to the motion and proportional to 
some power of the speed of the body relative to the fluid. The resistance to the motion 
of the bag of coins through the air is an example of this effect, as is the air resistance 
acting on moving vehicles. 


In Unit 4 two models for air resistance were investigated. They assumed the magnitude 
of the resistive force to be proportional respectively to speed, or to (speed)?. In 
mathematical modelling it is a good principle to start with the simplest available model, 
and not to complicate matters unless this initial model turns out to be inadequate. We 
shall therefore model the air resistance on the bag of coins as being proportional to the 
speed of the bag (which, since the air is substantially stationary, is also the speed of the 
bag relative to the air). Thus if the velocity of the bag is < then the damping force is 
assumed to have magnitude R = r|z|, where r is a positive constant, and direction 
opposing the motion. This is a much used model called linear damping. It has the 
virtue of simplicity and is mathematically easy to handle, which is such an advantage 
that the model is sometimes used even in situations where it may not be strictly 
appropriate. Any part of the frictional force on the bag of coins which is due to effects 
other than air resistance (such as internal friction in the elastic string) is a case in 
point. In such circumstances one would not expect the model to give results which are 
accurate in detail, but it is capable of showing the effects of friction in general terms. 


Linear damping 


We use a linear model for damping which stipulates a damping force of 


‘magnitude R = r|z|, where & is the velocity of the particle and r is a positive The SI units of the damping 
constant. The direction of the damping force is opposite to that of the constant r are Nm ‘s (force 
velocity. divided by speed), or 
equivalently kgs~’. 


1.3. The revised equation of motion 


We now investigate how the inclusion of linear damping alters the predictions derived 
from our model for the oscillating bag of coins, and the first step is to obtain the new 
equation of motion. As in Unit 7, the bag of coins is modelled as a particle of mass m 
and the elastic string as a perfect spring of stiffness k. Here, however, the particle is 
also acted upon by a damping force of magnitude R. Figure 3(a) shows the model 
without specifying any forces. Its purpose is to establish the origin for the displacement 
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x of the particle. The most convenient origin for our purposes is the rest position of the 
particle, that is, the position in which it will remain stationary when released. 


Figure 3 


Exercise 1 (revision) 


A particle of mass m is hung from a perfect spring of stiffness k and natural length lo, as You tackled this problem as 
shown in Figure 3(a). Find the length 1; of the spring when the particle is in equilibrium. Exercise 9 of Unit 7 


Solution on page 39 | Section 2. 


As well as choosing the origin at the equilibrium position, the direction of the x-axis is 
taken to be downwards, as indicated in Figure 3(a). This means that: 


(i) «x is positive when the particle is below its equilibrium position and negative when 
it is above; 


(ii) & is positive when the particle is moving downwards and negative when it is 
moving upwards; 


(ili) < is positive when the acceleration of the particle is directed downwards and 
negative when the acceleration is directed upwards; 


(iv) the x-component of any force acting on the particle is positive if the force is 
directed downwards and negative if the force is directed upwards. 


Figure 3(b) shows the forces which act on the particle when it is moving downwards 
and the spring is extended. In accordance with the convention used in Units 4 and 7, 
the arrows in Figure 3(b) indicate the directions of the forces and the quantities written 
beside the arrowheads are the force magnitudes. Thus T represents the magnitude of 
the spring force; the spring is in tension, and so it exerts an upward pull on the 
particle. The damping force of magnitude R also acts in the negative x-direction (since 
its direction is opposite to that of the velocity), and the gravitational force on the 
particle is downwards, as always. Hence, adding the z-components of these forces and 
using Newton’s second law, we have 


mg -T—-R=mi. (1) 


We next seek to replace the quantities T and R by expressions which involve the 
variables x and x. The result of Exercise 1 shows that when the particle is in 
equilibrium (where x = 0) the spring has extension mg/k. If the particle has position x 
below its equilibrium position then the spring will be further extended by the distance 
x, so that its total extension for this particle position is mg /k +a. The magnitude of 
the spring force will then be the product of its stiffness and extension, that is, 


T=k(“2 +2) =mg+ kz. 


Also, since @ is positive (with the particle moving downwards), the model of linear 
damping introduced in Subsection 1.2 gives 


A= riet= 1: 
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Substituting these expressions for T and R into Equation (1), we obtain 
mg —(mg+kax) —rz = mi, 
or miz+rz+kz=0. (2) 


This is the new equation of motion, including the effect of linear damping. Notice that 
all the terms on the left-hand side have the same sign, and that if damping were to be 
removed (by putting r equal to zero) then the remaining equation would be that of the 
simple harmonic model derived in Unit 7. 


Equation (2) was derived by considering a particle moving downwards and an extended 
spring. However, it applies also in other circumstances. Suppose, for example, that the 
spring is still extended but that the particle is moving upwards. Then the situation will 
be as shown in Figure 4 below. 


Figure 4 


The direction of the spring force in Figure 4 is upwards as before, and the expression 
derived earlier for TJ in terms of x will be unchanged. The velocity z is now negative, so 
the damping force is directed downwards (in the positive x-direction) and has 
magnitude 


R= rf) = Fz: 

Thus Equation (1) is replaced by 
mg—-I+R=mz, 

but when —rz is substituted for R and mg + ka for T, we have 
mz+re+kx =0, 


which is Equation (2) once more. This demonstrates an important feature of the linear 
model for the damping force: the same equation holds throughout the motion, and in 
deriving this equation we need consider only one particular configuration. Usually it is 
simplest to assume that the velocity is positive, or equivalently, that the particle is 
moving in the positive x-direction. Similarly, as with the simple harmonic model in 
Unit 7, there is no loss of generality in assuming that the spring is extended. The case 
when the spring is compressed leads to exactly the same equation of motion. 


Exercise 2 


Redraw Figure 3 for the case in which x measures the displacement of the particle from its 
equilibrium position and the z-axis is directed upwards. For the force diagram, assume that 
the spring is extended and that the particle is moving upwards. Show that the corresponding 
equation of motion for the particle is Equation (2) above. 


[Solution on page 39 | 
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The analysis above used the equilibrium position of the particle as the origin of 
position, for the sake of simplicity. It is, of course, possible to use other origins, and it 
is sometimes preferable to measure the displacement of the particle from the fixed end 
of the spring, as shown in Figure 5. You saw a similar alternative for the undamped 
perfect-spring model in Unit 7. 


Example 1 


A particle of mass m is hung from a fixed point by a perfect spring of stiffness k and 
natural length lo. The particle is subjected to a linear damping force with damping 
constant r. 


(i) Find the equation of motion of the particle when its displacement x is measured 
downwards from the fixed end of the spring. 


(ii) Hence find the length J, of the spring when the system is in equilibrium. Figure 5 


Solution 


(i) The system is illustrated in Figure 5. We assume that the spring is extended and 
the particle is moving downwards, so that the velocity « is positive. (Other cases rT 
lead to exactly the same equation of motion.) The forces acting on the particle, as 
shown in Figure 6, are: 


(a) the gravitational force, which is directed downwards and has magnitude mq; 
™m 


(b) the spring force, which is directed upwards (since the spring is extended) and 
has magnitude T = k(x — lg); 


(c) the damping force, which is directed upwards (since the particle is moving 
downwards) and has magnitude R = rz. 


Newton’s second law gives r mg 
mz=mg—-TIT—R Figure 6 
= mg — k(x —lIo) — rg, 
or mi+re+kr=mg-+ kilo, 
which is the required equation of motion. 


(ii) When the particle is in equilibrium, we have ¢ = 0 and % = 0. Substituting these 
into the equation of motion above produces 


kx = mg + klo 
mg 
= — ++ [o. 
or x i + lo 


Hence the equilibrium length of the spring is 11 =l9 +mg/k. O 


Example 1(i) shows that when the particle displacement is measured from the fixed end 
of the spring, the corresponding equation of motion is 


mi+réz+kx =mg-+ kilo. (3) 


Notice that the left-hand side of this equation is identical in form to that of 

Equation (2). However, Equation (2) is homogeneous whereas Equation (3) has a 
non-zero right-hand side. A particular solution of the inhomogeneous linear differential 
equation (3) is given by 


kip = mg + kilo, 
or Ln = + Io, 


which corresponds to the equilibrium position of the particle. This confirms that 
Equations (2) and (3) are alternative descriptions of the same physical behaviour. 


Indeed, if the equation of motion is known with respect to one choice of origin, then we 
can find the equation of motion with respect to any other such choice by a change of 
variable. For example, the equation of motion of our damped vibrating system when 
the displacement of the particle is measured from the fixed end of the spring is given by 
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Equation (3). The displacement of the particle from its equilibrium position, measured 
downwards, is (see Figure 7) 


X=2r2-lh=2- (to + =). 
k 
From this we have r = X +19 + mg/k, ¢ = X and # = X. Equation (3) can then be 
written in terms of X as 


mX +rX +k(X+lo+—*) = ig. + tile, 
or mX +rX+kX =0, 
which is just Equation (2) with the variable X in place of z. 


Exercise 3 
The equation of motion (Equation (2)) for the damped vibrating system considered above is 
mzi+rz+k2=0 


when the displacement is measured downwards from the equilibrium position. Let X be the 
upward displacement of the particle from floor level (see Figure 8). Assuming that the 
equilibrium position of the particle is a distance h above the floor, use a change of variable to 
derive the equation of motion in terms of X. 


[Solution on page 389 | 


1.4 The condition for oscillation 


When the bag of coins considered earlier in this section is displaced from its equilibrium 
position and released, it oscillates with a progressive decrease in amplitude. Imagine 
now what would happen if, instead of being suspended in air, the bag were suspended 
in treacle or in thick oil. Intuitively, you can see that its motion would be much more 
sluggish and that vibration might not occur at all. The fact that the system may or 
may not oscillate, depending on the amount of damping present, is reflected in the 
types of solution which are obtainable from the equation of motion. 


The general form of this equation for the model we are considering is 
mi+rz+kxr= kre, (4) 


where z is the position of the particle relative to some origin and z, is the displacement 
of the equilibrium position from that origin. A mechanical system whose equation of 
motion is of this form is called a damped linear harmonic oscillator, or damped 
harmonic oscillator for short. 


In Unit 6 you saw that the general solution of a differential equation such as 
Equation (4) is the sum of any particular solution and the complementary function, 
which is the general solution of the associated homogeneous equation 


mi+réz+kr=0. | (5) 


A particular solution of Equation (4) is x = x,, corresponding to the situation in which 
. the particle stays at rest in its equilibrium position. To find the condition under which 
the particle may oscillate about this point, it is necessary to solve Equation (5). Indeed 
if xe is taken to be zero, that is, if the origin for x is chosen to be at the equilibrium 
position of the particle, then Equation (5) is itself the equation of motion for the 
system. For simplicity we shall assume for the moment that this is the case, but note 
that the general solution of Equation (4) may be obtained directly from that of 
Equation (5) by addition of the term z,. 


The auxiliary equation corresponding to Equation (5) is 
mA? + rr +k =0, 
whose solutions are given by 


—r+V/r2 —4mk 


2m 


A= (6) 


If the solution for x in terms of time t is to be oscillatory then the roots \ must be 
complex, for which the condition is r? < 4mk. In other words, if the damping constant 


Figure 7 


Figure 8 


Here we are applying 
Procedure 1.1 of Unit 6. Note 
that this procedure leads to 
an oscillatory solution only in 
case (iii), when the auxiliary 
equation has complex roots. 


MST204 8.1 11 


r is small enough (for example, in air) then the system can vibrate, but if r is too large 
(for example, in treacle) then no oscillation occurs. 


If r? < 4mk then Equation (6) can be written as 


ae 


2m 


A= 


The complementary function is therefore 


V4mk — r? 
z(t) = Ae~"*/™) cos (ee + : 
m 


= Ae” cos(Nt + ¢), | (7) 


= — () = 
i, hi 


where 


(8) 


and A, ¢ are constants (with A non-negative and ¢ between —7z and 7) whose values 
may be chosen to fit the initial conditions of the system. Equations (7) and (8) give the 
general solution of Equation (5) provided that r? < 4mk. Equation (7) represents an 
oscillatory motion with angular frequency {2 and a diminishing amplitude Ae. FP’. 
The general shape of its graph is shown in Figure 9, which resembles in all essential 
respects the observed motion shown in Figure 2. Thus the inclusion of damping in the 
model for the motion of a bag of coins has enabled us to predict an important feature 
of the motion which did not show up in the simple harmonic model. The revised model 
is therefore an improvement on the first one. 


Figure 9 


The period of the damped motion is T = 27/Q, with 2 given by Equations (8). The 
corresponding period for undamped motion was found in Unit 7 Section 2 to be 27/wo, 
where wo = Wk /m is the angular frequency of simple harmonic motion. Not 


surprisingly, the expression (8) for 2 reduces to \/k/m when the damping constant r is 
equated to zero. 


Exercise 4 


Show that, for a given mass m and spring stiffness k, the period of linearly damped oscillations 
is greater than the period of undamped oscillations. 


[Solution on page 39 | 


The amplitude of the oscillations described by Equation (7) decreases steadily due to 
the exponential factor e~°’, and the greater the damping (that is, the larger the value 
of r), the faster this amplitude decreases. To verify this, we shall look at the ratio 


x(t +7) 
os. 3 aaa 


According to Procedure 1.1 
of Unit 6, the right-hand side 
of Equation (7) would have 
been written as 


e *(BcosNt + Csin Nt), 


where p, §} are as in 
Equations:(8) and B, C are 
arbitrary ‘constants. | 
However, we have written the 
term in brackets here in the 
alternative form given by - 
Equation (5) of Unit 7 
Section 2. 


In the remainder of this unit 
we use wo to denote the 
angular frequency of simple 
harmonic motion. The 
symbol w, which was 
employed for this purpose in 
Unit 7, will be given a 
different meaning in 

Section 2. 
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where T = 277/Q) is the period of damped oscillations. Using Equation (7), we have 
a(t +7) = Ae~°"*7) cos (Q(t + 7) + 6) 
= Ae +7) cos(Mt + 2x + 6). 
Now cos(Qt + 27 + ¢) = cos(Nt + ¢), and so 


—pe>T. 
x(t + T) = Aeelt+T) ee ga) oem Po (9) 


ce) ee ee 
As p > 0 and T > 0, we have 
ee 1 ame 2 ee ee. 


Equation (9) gives e~°” as the constant ratio of any two non-zero values of x one period 
apart; it applies, for example, to two consecutive maximum displacements on the same 
side of the rest position. As the damping constant r increases, it can be seen from 
Equations (8) that p increases while 2. decreases. Hence 7 increases and e °” decreases, 
confirming that increased damping produces a faster decrease in amplitude. The box 
below summarizes the results of this subsection. 


Linearly damped vibrations 

The equation of motion of a damped harmonic oscillator, namely 
mz+rae+kzx = 0, 

has an oscillating solution provided that 
r? < 4mk. 

The general solution in this case is 
a(t) = Ae~** cos(Qt + @), 


where 


es Gg. Ve 
am 2m’ a 2m 


and A, ¢ are arbitrary constants. All solutions of this form may be obtained 
with A non-negative and ¢@ between —7 and 7a. 


The angular frequency of these damped vibrations is 2, and their period is 
T = 27/Q. The ratio of two non-zero values of x one period apart is 


x(t +7) es 
x(t) 


We next return for a final look at the experimental results obtained from the suspended 
bag of coins, showing how measurements of the amplitude can be used to estimate a 
value for the damping constant r. 


With 20 coins in the bag, it was found that the maximum displacement of the bag from 
its equilibrium position decreased by a factor of about 10 over 8 complete cycles. 
Taking x to represent displacement from the equilibrium position, this means that 
x(t + 87) 
x(t) 
Now the left-hand side of this equation can be re-expressed as 


a(t+8r) x(t+77r) x(t +7) = (e~P7)8 = 87, 


a(t+%r) w(t+6r). |... x(t) 


ee 1) 1. 


Thus we have e~8°7 ~ 0.1, from which 


pt ~ —% log, 0.1 = 0.288. 
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Using Equations (8) and 7 = 277/Q) gives 


r 2n x 2m 
— X ——————_. ~ 0.288, 


2m /4mk — r? 

r? = 0.2887 
ee ae 
This equation may be solved to obtain a value for r provided that numerical values for 
the constants m and k are available. In Unit 7 Subsection 2.4 it was stated that each 
coin has a mass of 0.01 kg, so 20 coins have a mass of 0.2kg. It was shown also that the 
elastic string can be modelled by a perfect spring of stiffness 2.28Nm~!. Solution of 
the equation above with m = 0.2 and k = 2.28 results in the estimate 


r~0.062Nm—!s 


Or 


for the damping constant. 


As you saw above, the condition for oscillation is r? < 4mk. For the current case, this 
means that oscillation is possible provided that 


p< 402K 298 or r= 1.35. 


The estimate 0.062 for r is less than a of this upper bound, showing that quite a small 
damping force is enough to make the amplitude decrease fairly rapidly. The effect of 
small damping on the angular frequency and hence on the period is less marked, as you 
- are asked to demonstrate in the next exercise. 


Exercise 5 


Estimate the angular frequency of oscillations when there are 20 coins in the bag, as predicted 
by 


(i) the model for linearly damped motion, using the value of r derived above; 


(ii) the model for undamped motion developed in Unit 7. 


Exercise 6 


If as before m = 0.2, k = 2.28 and r = 0.062, find the values for A and ¢ in Equation (7) which 
correspond to the initial conditions x(0) = —0.05 and z(0) = 0. (Recall that A should be 
non-negative and that ¢ lies between —7z and 7.) 


Exercise 7 


Figure 10 shows a highly simplified model of a car. The particle of mass m represents the car 
body. This is supported by a perfect spring of stiffness k and a ‘damper’, which together 
represent the suspension. The contribution of the wheels is neglected and only the vertical 
motion is modelled. The damper contains a piston which moves through oil inside the 
cylinder. This causes a resistance to the motion of the piston which is approximately 
proportional to its velocity relative to the cylinder, the constant of proportionality being r. 


car body 


As part of a quick test of the suspension by a potential buyer when the car is stationary on a damper 


horizontal road, the body is depressed 0.1m below its rest position and then released. If 

m = 800, k = 5 x 10* and r = 9 x 10°, all in the appropriate SI units, establish whether the 
ensuing motion is oscillatory and, if so, find the distance of the car body below its rest position Figure 10 
after one complete cycle. 


[Solutions on page 39 | 


1.5 Weak, critical and strong damping 


In Subsection 1.3 we derived the equation of motion 
mz+rz+kxr=0 


for a mechanical system consisting of a particle of mass m moving under the influence 
of a perfect spring of stiffness k and linear damping with damping constant r, where 
the displacement x of the particle is measured from its equilibrium position. In 
Subsection 1.4 the general solution of this equation was found for the case in which the 
inequality r? < 4mk is satisfied. There are, of course, solutions also for the cases when 
r? = 4mk and r? > 4mk, and you are asked to derive these in the next exercise. 
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Exercise 8 

Find the general solution of the differential equation 
mz+rze+kax = 0, 

where m, r and k are positive constants, when 


(i) 1? <4mk, Part (i) of this exercise is 
isl ion 1.4. 
(ii) pik revision of Subsection 


(iii) r? > 4mk. 
[Solution on page 40 | 


The results of this exercise are summarized below. 


Weak, critical and strong damping 
The equation of motion 
miz+rz+kxr =0 
has the following solution. 
(i) Ifr? < 4mk (weak damping) then the general solution is 
a(t) = Ae~** cos(Mt + ¢), 


where 


rrr ea 


2m 


pe sees 
2m 


Here A and ¢ are constants; all solutions of this form may be obtained 
with A non-negative and ¢ between —7 and 7. 


If r? = 4mk (critical damping) then the general solution is 
gil = Be iC ©. 

where p = r/2m and B, C are arbitrary constants. 

If r? > 4mk (strong damping) then the general solution is 
sii Be Ste 

where B, C are arbitrary constants and 


pofela—4mk dil? 4k 


a= Fess NOESY eis a 


The case r? < 4mk, which is called weak damping, was discussed in the previous 
subsection. The general solution 


a(t) = Ae cos(Nt + ¢) 


is the product of an exponential decay and a sinusoidal oscillation of angular frequency 
Q. So the particle oscillates backwards and forwards about its equilibrium position with 
period 7 = 27/Q, much as it would when there is no damping, but with the amplitude 
of the oscillations decaying exponentially. The graphs of typical undamped and 
weakly-damped oscillations are shown in Figure 11 below. 


When r? > 4mk, which is called strong damping, the solution is 
a(t) = Be~*!* + Ce 2", 


where p; and p2 are positive constants. This solution is the sum of two decaying 
exponential terms. So from any initial conditions the particle returns to its equilibrium 
position without oscillation; the damping force dominates, and prevents any vibration. 
A graph for the strongly-damped case appears in Figure 11 below. 


In the case r? = 4mk, which is called critical damping, the solution is 


a(t) = Be? + Cte”, 


MST204 8.1 15 


where p is a positive constant. Both parts of this solution represent a decaying motion 
and so the motion again fails to be oscillatory, as shown in Figure 11. Given any initial 
displacement and velocity, the particle returns asymptotically to its equilibrium 
position just as for strong damping. As suggested by the figure, critical damping gives 
the fastest return to the equilibrium position. Hence many physical devices, such as 
recoil mechanisms for large artillery guns, mechanical suspensions, door-closing 
mechanisms and weighing scales, are designed to have critical or near-critical damping. 


strongly 
damped 


critically 


~ 
~ 
~ 
~, 
~ 
~ 
ign 
~ 


_—— 
- 
Soe ee 
_—— 


ae oe 
= damped 


undamped 


Figure 11 


Exercise 9 

A gun has a barrel, of mass 10° kg, which starts to recoil at a speed of 30ms~! when the gun 
is held horizontally and fired. Immediately after firing, the barrel starts to compress a spring 
which brings the barrel instantaneously to rest after it has recoiled a distance of 2m. In this 
phase of the barrel’s motion no damping is present. A linear damping mechanism is then 
engaged to ensure that the barrel returns to its initial position without overshooting the mark. 
(i) What is the spring stiffness? 

(ii) What damping constant is required for critical damping? 


[Solution on page 40 | 


Summary of Section 1 


1. Damping can be modelled by a damping force of magnitude r|z|, where r is a 
positive constant and z is the velocity of the particle, and direction opposed to the 
velocity of the particle. The SI units of the damping constant are Nm7's (or 
equivalently kgs‘). 


2. The equation of motion of a particle of mass m moving under the influence of a 
perfect spring of stiffness k and linear damping with damping constant r is 
mz+rz+kz =0, 


where z is the displacement of the particle from its equilibrium position. Such a 
system is called a damped harmonic oscillator. 


3. For weak damping (r? < 4mk) the general solution of the equation of motion is 
a(t) = Ae~™ cos (Nt + ¢), 


where 


r V4mk — r? 
p=— and Q = ———. 

2m 2m 
All solutions of this form may be obtained with the arbitrary constants A and ¢ 
restricted so that A is non-negative and ¢ lies between —7 and xz. This motion is 
oscillatory with period 7 = 27/Q. The amplitude of the oscillations decreases 


exponentially, so that 


Sit + Tet): ce P* 


16 MST204 8.2 


4. For critical damping (r* = 4mk) the general solution is 
a(t) = Be~* + Cte—™, 


where p = r/(2m) and B, C are arbitrary constants. This case gives the fastest 
return to the equilibrium position. 


5. For strong damping (r* > 4mk) the general solution is 
eG) = Be 74+ Ce 9. 
where B, C are arbitrary constants and 


— r+: fe? — dink 


pi ) Pa = 


2m 


2 Forced vibrations 


2.1 Introduction 


So far in this unit we have looked at the behaviour of a damped harmonic oscillator 

which moves freely after being started off from an initial position and with an initial 

velocity. We now go on to consider the motion of such a system when it is maintained 

in a state of vibration by some external physical agency. This ‘agency’ is modelled by a 

force changing continuously with time. The model of a car considered in Exercise 7 of 

Section 1 suggests one example: if the engine is started with the car stationary then it 

exerts an alternating force on its supports and hence on the car body, which vibrates 

for as long as the engine is running. As another example, think of a voltmeter which is 

used to measure a voltage varying with time. As the voltage increases and decreases, 

the pointer of the instrument will swing backwards and forwards. In both of these 

cases, the external force applied to the vibrating mass may vary with time in quite a 

complicated manner. We shall assume that the forcing agent is periodic and in that 

case, as will be shown later in this course (Unit 31), the force can be expressed as the 

sum of a number of sinusoidal terms. It follows that if we know the effect of a 

sinusoidal force we can, at least in principle, cope with any periodic variation by This ‘superposition of effects’ 
superimposing the effects of as many sinusoidal components as necessary. In this is dealt with by applying the 
section, therefore, we shall concentrate on the motion of the simple damped model superposition principle of 


se : : Unit 6 Subsection 2.4 to 
when it is subjected to a sinusoidal force. pubes eudiial Gauction 


involved. 


2.2 The equation of motion 


t & 
R 
£ 
Pcoswt 
mg 


Figure 1 
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Figure 1(a) shows the model system of Figure 3 of Section 1, where the origin for 
position x is chosen to be at the particle’s equilibrium position and the direction of the 
z-axis is taken to be downwards. This system is as before acted upon by the force due 
to gravity, a perfect spring and linear damping, but now there is in addition an applied 
sinusoidal force with x-component Pcoswt. The positive constant P is the amplitude 
(maximum magnitude) of the applied force, and w is its angular frequency. These four 
forces are shown in Figure 1(b) for the case in which the spring is extended, the particle 
is moving downwards and P cos wt is positive. However, the same equation of motion 
results also for the other possible configurations of the system. 


As before, the magnitudes of the spring and damping forces for the case shown are 
respectively 1’ = mg + kx and R= rz. Application of Newton’s second law to Figure 1 
therefore gives the equation of motion 


mg + Pcoswt — rz —(mg+kxr) = mi, 
or Mase + ke = PF coswt. (1) 


As you may recall from Unit 6, the general solution of this equation is the sum of two 
parts: the complementary function, which is the general solution of the associated 
homogeneous equation (with a zero right-hand side), and a particular solution of the 
original equation, which depends upon the form of this equation’s right-hand side. The 
complementary function for Equation (1) was found in Section 1; in this section we 
concentrate on the particular solution. 


When, as here, the equation models a physical system (with m, r and k positive), the 
amplitude of the complementary function decreases steadily and eventually becomes 
negligibly small. For this reason the complementary function is often called the 
transient part of the solution of Equation (1). On the other hand, the particular 
solution is persistent since, as you will see, it has a constant amplitude. It therefore 
remains in evidence long after the transient part of the solution has become negligible, 
which explains why the particular solution is known as the steady-state part of the 
solution. 


Since the right-hand side of Equation (1) is sinusoidal, the easiest way to find its 
steady-state solution is to use the phasor method. This is based upon Euler’s formula, 


0 


e’” = cos@ +7sin 6. 


Putting wt in place of @ gives 


et — coswt + isinuwt, 


so that cos wt is the real part of e’“*. It follows that Pcoswt = Re(Pe™*). A particular 
solution of Equation (1) is now obtained by trying a solution of the form 
x = Bcoswt + C'sinwt, or equivalently, 

x = Acos(wt + ¢). 


This is sinusoidal, of the same frequency as the external force but out of phase with it 
by an angle ¢. Application of Euler’s formula with wt + ¢ in place of 6 leads to 


2 = Re(Ae‘“*+?)) — Re( Ae**e™*). 


We next put Ae’? = Z, where Z is a complex constant (the phasor of the sinusoidal 
solution). The modulus of Z is A, the amplitude of the steady-state solution. The 
argument of Z is ¢, the phase of the steady-state solution. For x and its first two 
derivatives we then have 


= Re( Ze"): (2) 
& = Re(iwZe™*), 
B= Rett w Ze") = Re(—w* Ze). 

Substitution of these expressions into Equation (1) produces 


Re[(—mw? + irw + k)Ze™*] = Re( Pe). 


Note our notation for the 
three types of angular 
frequency occurring in this 
unit: (i) wo for simple 
harmonic motion (undamped 
and unforced); (ii) Q for 
damped harmonic motion 
(unforced); (iii) w for an 
applied force (and for the 
steady-state oscillation which 
it causes). 


See the box on page 14. 


This behaviour was | 
illustrated in Figure 11 of 
Section 1. 


See Procedure 2.2(a) in 
Unit 6. 


Unit 5 Subsection 4.1 
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Hence Equation (2) will represent a solution of Equation (1) if 


(—mw? + irw+k)Z =P, 


giving 
ae a ee F 
— —mw? +irwtk (k— mw?) +irw 
ag ae 
= oS P|(k — mw*) — irw] 


(k — mw?)? + r2w? © 


This is the phasor of the steady-state solution for x. Its relationship to P (the phasor 


of the applied force) is shown graphically in Figure 2. It may be expressed as Figure 2 
Z = u-— iv, where 
P(k — mw?) Prw 
u = — > ~— SO and ov = —_.__.. 
(k — mw?)? + r2w? (k — mw?)? + r2w? 


Hence the modulus and argument of Z are respectively 


P 
A=) est? = eee (3) 


(k — mw?)? + r2w? 


and o = Arg(Z) = —arccos (a) = — arccos (es | (4) 


k — mw?)? + r2w? 
The steady-state solution of Equation (1) is therefore 
x = Acos(wt + ¢), (5) 


where A and ¢ are given by Equations (3) and (4) respectively. This sinusoidal 
solution, with amplitude A, is the ‘output’ from the system corresponding to the 
‘input’ supplied by the forcing term with amplitude P. The angular frequency w of the 
output is the same as that of the input, but between input and output there is a time 
delay of —¢/w, as illustrated in Figure 3 below. 


Figure 3 


Exercise 1 
Derive Equations (3)—(5) without the use of phasors, by substituting the trial solution 
x = Bcoswt+ C'sinwt 


into Equation (1). (This is the alternative method of Procedure 2.2 in Unit 6.) 
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Exercise 2 


Consider once more the model car described in Exercise 7 of Section 1, for which m = 800, 
k=5x 10* and r = 9 x 10°, all in the appropriate SI units. While the car remains stationary 
on a horizontal road, its running engine produces a force on the car body which may be 
considered as sinusoidal, with amplitude 800 N and frequency 2.5 Hz. (Recall from Unit 7 
Subsection 2.3 that 1 Hz = 1 cycle per second.) 


(i) What is the corresponding value of the angular frequency w in rads !? 

(ii) What is the amplitude of the steady-state vibration experienced by the car body? 

(iii) What is the time lag between the maximum value of the engine force being exerted in the 
upwards direction and the car body next reaching its maximum displacement in that 
direction? 

Exercise 3 

The equation of motion of a mechanical system is 

2x + 52 + 142% = cost. 
Use the phasor method to find the steady-state oscillations of this system in the form 
x(t) = Acos(wt + @). 


[Solutions on page 40 | 


You have seen above how the phasor method can be used to find the steady-state 
vibrations of a model system when the external force has x-component P coswt. The 
method employs the relationship 


P coswt = Re(Pe**). 


A similar phasor approach can be adopted when the applied force is of the more general 
form P cos(wt + €), or of the equivalent form M coswt + N sinwt. In these cases we use 
the equations 


P cos(wt + e) = Re(Pe*e™*), 
M coswt + N sinwt = Re ((M — iN )e*) 
respectively. You are asked to consider the second of these possibilities in the exercise 
below. 
Exercise 4 
(i) Verify that 
Re ((M — iN)e'**) = Mcoswt+ Nsinuwt. 
(ii) Use the phasor method to find the steady-state solution of the equation of motion 
x+ 2x2 + 5x2 = 2cos 2t + 9sin 2¢, 
giving your answer in the form 
z(t) = Mcoswt + N sinwt. 
What are the amplitude and phase of this vibration? 
[Solution on page 41 | 


2.3 The damping ratio 


In the previous subsection we used the phasor method to show that the equation of 
motion of the forced and damped harmonic oscillator, 


mi+rz+kxz = Pcosut, (1) 
has the steady-state solution 

x = Acos(wt + ¢), | (5) 
where 

eat eae (6) 


(k — mw?)? + r2w? 


ae 2 
and d@ = — arccos ( | (7) 


(k — mw)?.+ r2w? 
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The applied force here is completely described by two parameters, its amplitude P and 
angular frequency w. For a particular set of values of m, r and k, the steady-state 
solution ‘x is a function of both P and w, and we wish now to examine the nature of 
this functional dependence. This amounts to investigating how the output parameters 
A and ¢ depend upon the input parameters P and w. As has already been noted, the 
angular frequency w of the output is identical to that of the input. 


Notice first that the phase angle ¢ is independent of the input amplitude P while, for a 
fixed value of the input frequency w, the output amplitude A is proportional to P. It is 
not so straightforward to see the effect upon A and ¢ of varying either the angular 
frequency w or the amount of damping present, but this is the principal aim of the 
current subsection. 


We proceed by reorganizing the information contained in Equations (6) and (7), so that 
the quantities m, r and k are replaced by combinations which permit a clearer view of 
how A and ¢ react to changes in w. In order to provide some motivation for the choice 
of these ‘combinations’ of m, r and k, we return first to Equation (1). Dividing this 
equation by m gives 


a P 
£+—x£+ —zxr = — cosut. (8) 
m m m 


The coefficient of x is now k/m, a quantity encountered previously. In Unit 7 
Subsection 2.3 you saw that the angular frequency of the unforced and undamped 
perfect-spring system was wo = ,/k/m, so we may replace k/m by w?. In this context 
we shall refer to wo as the undamped angular frequency. 


The next step is to rewrite the damping term. You saw in Section 1 thatthe 
complementary function of the equation of motion is oscillatory if and only if the 
damping constant r is less than an upper limit given by r? = 4mk, or r = 2V/mk. This 
leads us to define the damping ratio a of the system by 

r 


= 2/mk 


The damping constant r provides an absolute value for the damping force per unit 
speed exerted on the particle. On the other hand, the damping ratio a gives a measure 
of the importance of damping relative to the mass m and spring stiffness k of the 
system. Note that a is a dimensionless quantity, since the dimensions of r (Nm~'s, or 


kgs~') are identical to those of Vmk. 


From the.way in which a is defined, and from the results of Subsection 1.5, it follows 
that a system moving in accordance with the complementary function of Equation (8) is 


(i) weakly damped (oscillatory) if a < 1, 
(ii) critically damped if a = 1, 
(iii) strongly damped if a > 1. 


Figure 4 shows graphs of the complementary function for a particular system for four 
different values of a; the initial conditions are the same for all of them. This is 
essentially a repetition of Figure 11 in Section 1. 


Figure 4 


As mentioned before, we use 
wo in this unit to denote the 
natural (unforced, 
undamped) angular 
frequency, rather than w as 
in Unit 7, since currently w 
represents the angular 
frequency of the forcing term. 
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From the equations 


ck 
m Ay Nik 


we have 


r. QJow/imk ; 
= SS = 204] — = Zerit, 
m m m 


so that Equation (8) can now be written as 
: 2 re 
E+ 2awopkt + wor = — coswt. (10) 
m 


Exercise 5 

Determine the phasor for the steady-state solution of Equation (10). Hence find expressions 
for the amplitude and phase angle of this motion. 

Exercise 6 


A linearly damped perfect-spring system with sinusoidal applied force has the equation of 
motion 


47 + 9x + 100x = 9 cos 6t. 
(i) Find the values of the undamped angular frequency and the damping ratio. 
(ii) Determine the amplitude, phase and angular frequency of the steady-state solution. 


[Solutions on page 41 | 


In Exercise 5 you showed that the steady-state solution of Equation (10) is 
x = Acos(wt + ¢), where 


P/m 


he eS Sere eee 11 
(we — w?)? + 4a2wew? ay) 
2 8 
and @ = — arccos eg (12) 
(wg — w?)? + 40?2wew? 


These expressions could also have been obtained directly from Equations (6) and (7), 
using Equations (9) to replace k by mw and r by 2amwp. 


The forced and damped harmonic oscillator 
The equation of motion 
mi+ré+kxr = Pcoswt 


of the forced and damped harmonic oscillator can be rewritten as 
: ; D - 
L + 2awpt + wor = — coswt, 
m 


in terms of the undamped angular frequency 


k 
Wo = —— 
m 


and the damping ratio 


The steady-state solution z = Acos(wt + ¢) of this system has amplitude 
P/m 
/ (we — w?)? + 4a2wew? 


A= 


and phase 


2 
We — w 
o = —arccos | —————2¥ _________ ]} . 
J (we — w?)? + 4a2W2w? 


22 MST204 8.2 


The best way to show the significance of Equations (11) and (12) is to draw graphs. In 
the case of Equation (12), what we want to demonstrate is the effect on the value of ¢ 
of varying w and a. But the equation also contains wo, so one possibility would be to 
plot, for each of a range of values of wo, a set of curves of ¢ against w, each separate 
curve corresponding to a particular value of a. This might mean drawing a large 
number of sets of curves, which would require much calculation and plotting. In order 
to simplify the procedure, we note that Equation (12) can be rewritten as 


1 — (w/w)? 
(1 — (w/wo)?)? + 40? (w/w)? 


which shows that ¢ can be regarded as a function of the two variables a and w/wo. We 
therefore plot graphs of ¢ against w/we for particular values of a, as shown in 

Figure 5(b) below. This enables all of the required information to be presented by 
drawing just one set of curves. 


@ = — arccos ( 


The same sort of approach can be applied to Equation (11), except that the task of 
representing graphically the variation of A with w and a is further complicated here by 
the presence of m and P as well as wo. Noting that 


oa ees 2P 
aoe ae eS 


? 


and putting M = Ak/P, Equation (11) can be rewritten as 
1 


= /(1 — (w/wo)?)? + 4a?(w/wo)? 


where, as before, the right-hand side is a function of the two variables a and w/wp. The 
resulting graphs of M against w/w for various values of a are shown in Figure 5(a) 
below. 


(13) 


Notice that M, like a and w/w, is a dimensionless quantity. For a given spring stiffness 
k;, the ratio of the output displacement amplitude A to the input force amplitude P can 
be found from the value of M. 


(b) -m 


Figure 5 
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Exercise 7 


Calculate the value of M which corresponds to the steady-state solution of the differential 
equation given in Exercise 6 above. 


[Solution on page 42 | 


2.4 Resonance 


The most noticeable feature of the graph of M against w/wo in Figure 5(a) is that for 
low values of a it shows a definite peak near w/w o = 1. This means that, when the 
angular frequency w of the sinusoidal external force is close to the undamped angular 
frequency Wo of the spring—mass system, the steady-state amplitude of the displacement 
of the particle becomes large compared with its values for the same force amplitude 
over most of the frequency range. The motion of the system at these relatively large 
amplitudes is known as resonance. As a increases, the resonance peak becomes less 
pronounced and moves towards lower values of w/w. When a = 1/\/2, the peak is at 
w/wo = 0 (as you are asked to show in Exercise 10 below). For a > 1/,\/2 the curve has 
negative slope for all positive values of w/wo. 


Resonance is found in many physical systems and may be considered desirable, as in 
the case of a radio receiver which is tuned to an input signal of a particular carrier 
frequency, so that inputs at other frequencies are neglected. However, there are other 
situations where resonance may be most undesirable. In the case of a car, large vertical 
movements of the car body will certainly be uncomfortable and may also be dangerous 
and destructive. 


Figure 5(a) may be used to represent the response of an instrument which measures 
and records a physical quantity such as, for example, electrical voltage. The movement 
of the instrument’s pointer (and attached parts of the mechanism) can be modelled by 
the movement of a particle, and the voltage causes, and is proportional to, the external 
applied force. If the voltages to be measured are periodic and have a range of different 
frequencies, then it is preferable to avoid re-calibrating the instrument for every 
different frequency; indeed, it might be the case that the voltage consists of several 
sinusoidal functions of different frequencies superimposed. Hence the same value of M 
is required for as wide a frequency range as possible. Of the curves drawn in 

Figure 5(a), that for a = 0.6 meets most closely the condition that M should be 
approximately constant over the greatest possible frequency range. The value usually 
chosen in practice is a = 0.65. In order to achieve this damping ratio, the instrument 
movements often incorporate deliberately added linear, or nearly linear, damping of a 
mechanical or electrical nature. | 


In order for the output to be a faithful representation of the input we need to have the 
same value of M for all the frequency components, as already explained, so that each 
component of the input makes the appropriate contribution to the output. In addition, 
we need to ensure that the time shift between different components (the interval 
between the instants at which different components reach their maxima, say) is not 
affected by passing through the instrument. If it were, the output would be a distorted 
version of the input. So, to avoid distortion, one wants the time lag introduced by the 
instrument to be the same for each component. This time lag equals —¢/w (see 

Figure 3, page 18), so that we want the phase ¢ to be directly proportional to the input 
frequency w. As it happens, for a = 0.65 this is almost true over a frequency range 
similar to that for which M is approximately constant. Figure 5(b) illustrates this. 


Figure 5(b) also shows that, for all values of a and w/wo, the phase ¢ lies between —7 
and 0, and that ¢ = —7/2 when w = uw. 


Exercise 8 


A voltmeter is modelled by a forced and damped harmonic oscillator with damping ratio 

a = 0.6 and undamped angular frequency wo = 20rads~’. Estimate from Figure 5(a) the 
frequency range over which the instrumert will reproduce the amplitude of a sinusoidal input 
voltage with a relative error of less than 10%. 


Exercise 9 


Use the graphs in Figure 5(a) to check your answers to parts (ii) and (iii) of Exercise 2 
(page 19). 
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Exercise 10 


Write the right-hand side of Equation (13) in terms of a and p = w/wo. Using this expression 
show that, for a given value of a, the maximum value of M occurs 


(i) when p= V1 — 2a?, if a < 1/V2; 
(ii) when p=0, if a > 1/V2. 
[Solutions on page 42 | 


Summary of Section 2 


- 


The equation of motion of a damped harmonic oscillator subject to an applied 
sinusoidal force of amplitude P and angular frequency w can be written as 


mz+kz+raz= P cosut. 


The steady-state forced vibrations of this system can be found by using the 
phasor method. 


This equation of motion can be rewritten in the form 
: 2 P 
£ + 2awypt + wor = — cosut, 
m 


where the undamped angular frequency w is defined (as in Unit 7) by 


k 
| ae ee. 
m 


and the damping ratio a is defined by 


es 2/mk- 


The steady-state vibration of the forced and damped harmonic oscillator is given 
by 


x = Acos(wt + ¢), 


where 
tes 
A= 
(we — w?)? + 4a?w2u? 
and go = — arccos Bieta RRA 5) 
(we — w?)? + 4a2wew? 


For a < 1/\/2, the amplitude of the forced vibrations exhibit a maximum at a 
particular value of the forcing angular frequency w. This effect is known as 
resonance. 
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3 The perfect dashpot 


3.1 Introducing the dashpot 


In Section 1 we augmented the Unit 7 model of a particle-spring system by the 
inclusion of linear damping, and in Section 2 an external force acting upon the particle 
was added. This section shows how the linear model for damping can be applied to 
configurations of particle, spring and damping which differ from those encountered 
previously. These configurations arise from situations in which an external force acts 
upon a point of the system other than the particle itself. 


For example, consider the vertical oscillations experienced by a motor car as a result of 
its motion along an undulating road. This may be modelled by considering the car to 
be at rest with the road surface moving beneath it, the effect of the undulations being 
represented by a sinusoidal displacement of the wheels about the mean road level. In 
predicting the vibrations of the car body under these circumstances we need to include 
a linear damping force which is proportional to the relative velocity between the body 
and wheels, rather than a force proportional to the absolute vertical velocity as in the 
earlier model for a stationary car. Also, we now have to incorporate an external force 
described by the vertical displacement of the wheels rather than a force acting directly 
upon the car body. 


The main purpose of this section is to enable you to write down the correct equation of 

motion for the particle when the system consists of any configuration involving a 

perfect spring, linear damping and an applied force. As a first step, Figure 1 shows a 

diagrammatic means of indicating linear damping. It represents a piston inside a 

cylinder, or dashpot. One form of dashpot is used as the shock absorber on a car (see — |\— 
Exercise 7 in Section 1), where the cylinder is full of oil and the relative motion of the 

piston forces the oil to flow through the annular space between piston and cylinder. 

This produces a resistance to the relative motion. For modelling purposes we assume Figure 1 
that damping is produced by a perfect dashpot, for which the resistance is 

proportional to the relative velocity of the dashpot’s ends (or equivalently, to the rate of 

change of the dashpot’s length). The constant of proportionality is denoted as before 

by r, but is now referred to as the dashpot constant rather than the damping 

constant. If the instantaneous length of the perfect dashpot is / then the magnitude of 

the damping force produced is r{/|. 


r 


The occurrence of a perfect dashpot on a diagram indicates the presence of linear 
damping, irrespective of whether the damping being modelled is produced by an actual 
dashpot or by some other physical means. In this respect the role of the perfect 
dashpot is analogous to that of the perfect spring, which indicates the presence of a 
force depending linearly on the relative displacement of two end-points, whether or not 
the real force being modelled is applied by the agency of a spring. 


Figure 2 shows how we could draw a perfect dashpot alongside a perfect spring to 
represent the model of Section 1 for the suspended bag of coins. The damping effects 
here are due not to a physical dashpot but to air resistance and internal friction in the 
rubber string. Since the point of suspension is fixed, the damping force magnitude r|I| 
in this instance is equal to r|z|, where x is the displacement of the particle from any 
fixed point. However, this will not always be the case. Figure 2 


Figure 3 shows a perfect dashpot with both ends in motion. The positions of these ends 
with respect to some fixed point are denoted by zx and y, so that their velocities are 
respectively z and y. Since the dashpot length is 1 = x — y, the magnitude of the 
damping force exerted by the dashpot is r|i |=r|z —y|. Now if z > y, that is, if the 
dashpot is lengthening, then the force exerted on whatever is connected to the dashpot 


at either end is directed towards the dashpot. (Imagine that you are pulling the two => —- 
ends of the dashpot apart with your hands. The dashpot will resist by exerting forces | | 
trying to stop your hands moving apart. The force on your right hand is to the left, and : z ! 
the force on your left hand is to the right.) By the same token, if « < y, that is, if the | y 


dashpot is shortening, then the forces on objects attached to the ends will be directed 
away from the dashpot, again resisting any change to the length of the dashpot. Figure 3 
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Example 1 


Figure 4 


A particle is attached to one end of a perfect dashpot, as shown in Figure 4. The 
displacements of the particle and of the other end of the dashpot from a fixed point O 
are denoted respectively by x and y, where x > y. Find the x-component of the 
dashpot force acting on the particle when 


(i) the length of the dashpot is increasing} 
(ii) the length of the dashpot is decreasing. 


Solution 
The length of the dashpot is / = x — y, and so i=ar- Y. 


(i) When the length of the dashpot is increasing, l is positive and so |/| =1 = « — y. 
The magnitude of the dashpot force is therefore R = r|l| = r(« — y). In this case 
the dashpot force is directed towards the dashpot, as shown in Figure 5, trying to 
resist the increase in the length of the dashpot. Hence the x-component of the 
dashpot force is 


—-R=-r(z—-¥y). 


(ii) When the length of the dashpot is decreasing, | is negative and so 
\i| = -l = —(& — y). Hence the magnitude of the dashpot force is 
R=r'\l| = —r(«#—y). Again, this force tries to resist the change in the length of 
the dashpot, and so is directed away from the dashpot, as indicated in Figure 6. 
Therefore the x-component of the dashpot force is 


R=-r(¢z-y). O 


This example demonstrates that the x-component of the dashpot force takes the same 
form whether the length of the dashpot is increasing or decreasing. Consequently, the 
equation of motion of a particle attached to a dashpot will also have the same form in 
these two cases. We shall therefore usually assume, when deriving an equation of 
motion, that the dashpot length is increasing (so that the dashpot force is directed 
towards the centre of the dashpot). However, the resulting equation of motion will be 
valid also for the case when the length of the dashpot is decreasing. 


The perfect dashpot 


The perfect dashpot is a diagrammatic representation of a type of force, 
rather than of the mechanism which produces it. It represents a resistive force 
proportional to the rate of change of the dashpot’s length /. In other words, 
the force has magnitude R = r|l|, where r is a positive constant known as the 
dashpot constant. 


If the dashpot’s length is increasing (so that J is positive), the dashpot force is 
directed towards the centre of the dashpot. If the length of the dashpot is 
decreasing (so that / is negative), the force is directed away from the centre of 
the dashpot. In either case the dashpot force resists any change in the length 
of the dashpot. 


The equation of motion for a particle attached to a perfect dashpot can be 
derived by assuming that the dashpot is increasing in length. 


Figure 5 


Figure 6 
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Example 2 


Figure 7 


A particle of mass m is attached to a fixed point O by a horizontal perfect spring of - 
natural length /p and stiffness k. The particle is also attached by a horizontal dashpot 
of dashpot constant r to a point whose displacement with respect to the fixed point O 
is y, as Shown in Figure 7. Find the equation of motion for the particle, assuming that 
it is otherwise free to move in a horizontal direction. 


Solution 


In deriving the equation of motion we shall assume that the spring is extended and that 
the length of the dashpot is increasing, so that the forces act in the directions indicated 
in Figure 7. We shall measure the displacement x of the particle from the fixed point 
O. The extension of the spring is x — Ip, so the magnitude of the spring force is 


T = k(x — lp). 


The length of the dashpot is / = x — y, and so i=a- y. Now I is positive, since the 
length of the dashpot is increasing. Hence the magnitude of the dashpot force is 


R=rl=r(é—¥y). 


By Newton’s second law, the equation of motion for the particle is 


mz=-T-—R 
= —k(x —I9) —r(z — 9), 
or mi+reze+kx=kigp+ry. O 


Exercise 1 


Find the equation of motion for the particle in Figure 2, taking xr to be the downward 
displacement of the particle from the fixed ends of the spring and dashpot, and /p to be the 
natural length of the spring. Compare your answer with Equation (3) in Section 1 (page 9). 


Exercise 2 


r | 


m 
Pcoswt 


Figure 8 


Find the equation of motion for the particle in Figure 8, which is attached to a perfect dashpot 
and perfect spring (of natural length Jo) as well as being subjected to an external force whose 
x-component is P cos wt. 


[Solutions on page 42 | 


Consider once more the equation of motion derived in Example 2 for the situation 
shown in Figure 7. Note that, as it stands, this equation of motion cannot be solved to 
obtain an expression for the function x(t), since the ry term on the right-hand side is an 
unknown function of t. This situation can be altered by specifying the position function 
y(t) of the left-hand end of the dashpot. This left-hand end is referred to as the forcing 
point of the system, since it is the point at which an external force is applied. 
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There are essentially three different configurations for a single-particle system, 
connected to one perfect spring and one perfect dashpot, which has an external force 
applied at a point other than the particle: 


(i) the external force may be applied to the other end of the spring, with the other 
end of the dashpot held fixed; 


(ii) the external force may be applied to the other end of the dashpot, with the other 
end of the spring held fixed; 


(iii) the other ends of the spring and dashpot may be joined together, and the external 
force applied to this joint end. 


Example 2 looked at an instance of case (ii). We next consider an occurrence of 

case (iii), including specification of the forcing point motion. Provided that this motion 
is sinusoidal, the mathematics required to solve the equation of motion for the particle 
is much the same as that employed earlier. 


Figure 9 is a diagram of the highly simplified model of a car which was mentioned at 
the beginning of this section. This model is obtained by neglecting the mass of the 
wheels and concentrating the wheel suspensions into one pair of perfect spring and 
perfect dashpot elements. The spring and the dashpot are assumed to have the same 
length at all times, with their lower ends being at the same horizontal level. The ‘input’ 
y(t) is a sinusoidal displacement of the common lower end of the spring and dashpot. 
This input represents the undulations of the road, which the car is traversing at a 
certain speed. The ‘output’ x(t) is the vertical displacement of the particle of mass m 
which represents the car body. Both x and y are taken to be positive upwards and 
measured from a fixed horizontal level, which is the mean height of the road surface. 


Assuming that the spring is extended and that the dashpot is increasing in length, the 
directions of the forces on the car body are as shown in Figure 9. The extension of the 
spring is x — y — Io, and the rate of increase of the dashpot’s length is « — y. So the 
magnitudes of the spring and dashpot forces are respectively 


T=k«2«-y-—lob) and R=r(ec—y). 


The gravitational force, of magnitude mg, acts vertically downwards. The equation of 
motion of the car body is therefore 


mz=-—-mg—-—TI—-R 
—mg —k(a —y—Ip) —r(a#—-y), 


or mi+rae+kae =klo -—mg+ky+ry. (1) 


If the undulations of the road are sinusoidal about its mean level and the car travels at 
constant speed, then the forcing point motion will also be sinusoidal. If these forcing 
point oscillations have amplitude Ap and angular frequency w then y = Ao coswt, giving 


mi +rxz+kx = klo —mg + Ao(kcoswt — rw sinwt). (2) 


This equation for forced and damped motion is similar to those which you met in 
Section 2. Its steady-state solution has the form 


x(t) = 1p — mg/k + Acos(wt + 9), 


where lg — mg/k is the height at which the particle could remain in equilibrium if no 
external force were applied. This constant part of the solution is obtained by finding a 
particular solution of the equation 


mi+rze+kz = kilo — mg. 


The sinusoidal part of the solution is derived by finding a particular solution of the 
equation 


mi+rz+kx = Ao(kcoswt — rwsinwt). (3) 


Exercise 3 


Derive an expression for the phasor of the steady-state solution of Equation (3), and hence 
derive expressions for the amplitude A and phase angle ¢ of the forced vibrations described by 
Equation (2). 


[Solution on page 42 | 


Figure 9 
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Under certain circumstances the position x of the particle relative to some fixed point 
may be of less interest than its position z relative to the forcing point. Since z = x — y, 
it is a simple matter to obtain the appropriate differential equation for z directly from 
that for x. For example, the height of the car body above the road surface in the model 
of Figure 9 is z = x — y. On replacing x by z+ y in Equation (1), we obtain 


mz+rz+kz=klo — mg — mj. 


Given that y = Ao coswt, this equation for z(t) can be solved in the same manner as 
was adopted for Equation (2). 


3.2 Deriving equations of motion (Audio-tape Subsection) 


The audio-tape provides further practice in setting up equations of motion, analysing 
two model spring-dashpot systems which are slightly more complicated than those 
considered previously. 


Start the audio-tape when you are ready. 


Note that during the tape the 
symbol v (Greek ‘nu’) is used 
in place of w for the angular 

frequency of the forcing term. 
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> First System 
particle 


one of mass m 


sinusoidal 


pected k pom external 
, = , force 
Vy l, 
Horizontal motion only 
constants ry & fy Se 


Equation of motion ? ~~ What forces act ? <- Variables ° 


€ System at Rest 


Take x as particle position. 


rest position of particle: 


direction 
io ae 
x 


: 
| | 
| , | 
€&> With External Force 


natural 
length of 


Spring 


Assume 


Bait | What are 
ex>0 the forces... 
ex>0 oes ¢ magnitude ? 
ePcos vt >0 e direction? 


4 The Force Diagram 
0 


x 
kx <—___— | ate | COS YT external 
mM force 


ix<— <——_ /,Xx 
dashpot 
forces 


Equation of motion ? 
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5 The Equation of Motion 


x-components of forces are 


Newton's 2nd law: 
MX = xX—component 
of total force 


e external: Pcos vt 
6 spring: ~kx 
e dashpots: teat, x 


mx=Pcos vt -kx -r,x-x 


Re-arrange: mx +(r,+r,) x +kx = Pcos vt 


OF | er kx =P cos yt}. ee 


6 Turned through 90° 


does If 
differ ? 


pm ; 
mM 
O 
| Oo 
gravity 
P cos vt 


c Second System — at Rest 


two springs, 
one dashpot 


© In Motion 


Assume x >0. 
ESO, yx. 


rest posifion of 
particle: origin for x 


Forces balance: 
mg mg=T=k(L,- fo) 
eal a fo + mg/k. 


Frame 5 equation applies again. 


r2 


Springs are at 
their natural lengths 


What are 
the force 
A\ magnitudes ? 


© spring 1 — © spring 2 © dashpot 


extension = 4 extension = ee length change rate = a 
tension = eo tension = | force magnitude = ee 
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10) The Force Magnitudes 


extension = 7 extension =| Y-X 
© spring < = a © spring 2< | sEGED 


tension fension = 


eee <— length change rate = Diagram of 
e dashpo 
i force magnitude = forces ? 


® The Force Diagram 


———_—___f— __ > 
0 x 
dashpot force 6 nen 
ae 
feces ce raeremnrnremnrereomenmmarmt 


Ky x 


Spring forces 


Equation of motion ? 


ol ig la all 


® The Equation of Motion 


x-components of forces are Nontons 2nd law: 


e dashpot: -rx mx = x-component 
€ spring ie =Kaz . of total force 
® Spring Zs +k, (Ym) 


mx =~1x-k,x +ky (y-x) 


Re-arrange: mx +rx +(k,+kz,)x=k2y 


, eka te 
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The remaining exercise in this section again asks you to derive an equation of motion, 
and illustrates how the response graphs of Figure 5 in Section 2 (page 22) can be 
adapted for use in different circumstances. 


Exercise 4 


Figure 10 shows a model of a solid body being moved vertically through water. The perfect Z 
dashpot represents the resistance of the water, and the perfect spring represents the resilience 

of a towing rope. (Note that a perfect spring can reasonably be used as a model for a rope 

only if the rope is in tension, since a rope, unlike a spring, cannot sustain a compressive force. ) 


(i) Initially the top end of the spring is held at a fixed position. If the natural length of the r 
spring is /o, what is its length when the particle is in equilibrium? 


oo 


(ii) The top end of the spring is now moved along a vertical straight line, with displacement y 
measured upwards from its previous fixed position. The corresponding displacement x of 
the particle is measured upwards from its equilibrium position. By making suitable 
assumptions, draw the diagram of forces acting on the particle, and derive the equation 


of motion for the particle. 
a : Figure 10 
(iii) Suppose that k = 200 Nm7! and m = 4kg, and that the damping ratio is a = 0.2. 


Suppose further that the top end of the spring is given a sinusoidal displacement about 
its original fixed position, of amplitude 0.05 m and frequency 2Hz. By comparing the 
equation of motion derived in part (ii) with Equation (1) of Section 2 (page 17), use 
Figure 5(a) in Section 2 (page 22) to estimate the amplitude of the particle’s steady-state 
oscillations. 


[Solution on page 43 | 


Summary of Section 3 


1. A perfect dashpot represents a resistive force which is proportional to the 
relative velocity between two components of a mechanical system. 


2. The dashpot force has magnitude R = rl |, where / is the length of the dashpot 
and r is a positive constant known as the dashpot constant. If the dashpot’s 
length is increasing, then the dashpot force is directed towards the centre of the 
dashpot; if the length is decreasing, then the force is directed away from the centre 
of the dashpot. 


3. The equation of motion of a particle attached to a perfect dashpot can be derived 
by assuming that the dashpot’s length is increasing. This equation of motion 
applies also to the case where the dashpot is contracting. 


4 Off the record: resonance and 
damping (Television Section) 


Now watch the television programme. 


Read the following after viewing the programme. ee: 
ce) 


Allan Solomon introduced the programme by saying that to understand why the 

tone-arm of a record player can misbehave, we need to understand damping and forcing TV8 
in simpler mechanical systems. Thereafter the programme divided naturally into three 

parts, each of which considered the vibrations of a particular apparatus: 


Part 1: the unforced, damped vibrations of an apparatus modelled by Figure 1; 
Part 2: the forced vibrations of an apparatus modelled by Figure 2: 
Part 3: the forced vertical vibrations of a tone-arm, modelled by Figure 3. 


Part 1 of the programme showed how a simple mass—spring—dashpot system 
behaved, and introduced the model for it shown in Figure 1. David Broadhurst found 
that a dashpot containing air made little difference to the oscillations observed without 
a dashpot: the motion still resembled simple harmonic motion (see Figure 4(a)). With 
water in the dashpot the damping was more noticeable, though the motion was still 
oscillatory (Figure 4(b)). Strong damping occurred with oil in the dashpot, when there 
were no oscillations but just a simple decay to the equilibrium position (Figure A(c)). Figure 1 
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(a) air 


m 


(b) water 


VES ELLE ES 
i Figure 3 
Figure 4 Output traces from a simple mass—spring—dashpot system. 
It was then shown that the equation of motion for the model in Figure 1 was This is Equation (3) of 
6 Section 1, which you derived 
mi +rz+kx = mg + kilo, once more in Exercise 1 of 
Section 3. 


where /p is the natural length of the spring. This equation was solved, and graphs of 
the predicted motion were shown for three values of the damping ratio a = r/(2Vmk) 
(see Figure 5, in which xo = Ip + mg/k is the equilibrium position of the particle). For 
a < 1 damped oscillations were forecast, while for a > 1 an absence of oscillations was 
predicted. Critical damping, with a = 1, gave the fastest decay to the static position. 
The observed results were clearly in qualitative agreement with the model’s predictions. 
One quantitative test is the subject of the exercise below. 


Exercise 1 


In Subsection 1.4 it was shown that, for the model of Figure 1, the ratio of particle 
displacements from the equilibrium position at times one period 7 apart is 


elt) 2 ae Gf 2: 3 2 Be 7 oot 
x(t) Figure 5 
where p = r/(2m). It follows from this that fig ci ages ae es aes 
model system shown in 
t 
Bee) ane for 152). 5. — 
x(t) 
Thus if Ao, Ai, A2,... are the heights above the mean of successive peaks on the graph of z(t), 
and Aj, A,, A2,... are the depths below the mean of successive troughs (with Aj between Ao 


and Aj, say), then 
Ag Aee + ang AS = Ane 


from which We add the two previous 
eee 1\_—npr equations here because it is 
An Ag Ss oe easier to measure A, + A’, 
= - ( Ao + Ao ) aoe from the trace than it is to 
©\ An + At measure A, or Aj, 
individually. 


Using the data provided below for the water-filled dashpot, test this prediction of the model 
by drawing an appropriate graph. 


ie case 
roca + aay [ae [aa Pas [ow [Or 


[Solution on page 44 | 
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Part 2 of the programme dealt with the forced oscillations of an apparatus 
modelled by Figure 2. The design of this apparatus is indicated in Figure 6. 


3 ae hinged y 
= — pens ; inge 3 arms y 
output = 


Figure 6 The apparatus used in Part 2 of the television programme. 


The motor-driven cam made the bottom of the spring move up and down sinusoidally 
with a fixed amplitude and at a forcing frequency that could be varied. This was the 
input oscillation, 


y = yo + Ao cos wt, 
recorded by the lower of the two traces. The steady-state output oscillation, 
© = 2% + Acos(wt + ¢), 


was recorded by the upper trace. In addition to the spring connecting the two arms of 
the apparatus, there was a dashpot connecting the upper arm with a point at a fixed 
level. The geometry is such that the effective damping constant can be increased by 
sliding the dashpot to the left. 


Initially the output was studied for three different forcing frequencies, with the same 
amount of damping. The first frequency gave an amplitude ratio A/Ao greater than 
unity. At a higher frequency the amplitude ratio was even greater, but at a higher 
frequency still it decreased. Thus the system exhibited resonance, since the amplitude 
ratio had a maximum greater than unity (see Figure 7). 


output | 2A 


inpek PPPS LSS SP ALLO POOP OPT PATS FPASP PAE = TRA 


Figure 7 Input and output traces for the apparaes of Figure 6 with light damping, 
near the natural frequency. 


The apparatus was modelled by the system of Figure 2, leading to the equation of 
motion 
 mé+ré+ka = klp — mg + ky, (1) 


where /g is the natural length of the spring. With y = yo + Ap coswt, Equation (1) 
becomes 


mz +rz+kz = k(lo + yo — mg/k) + kAp cos wt. (2) 
This has steady-state solution 
L= 20+ Acos(wt.t 9), 


where rp = Io + yo — mg/k and 
A = ‘ 1 (3) 
— Ao (1 — (w/wo)?)? + (2aw/wo)? 
Graphs of the amplitude ratio from Equation (3) were provided for values of the 
damping ratio a between 0.1 and 1.0 (see Figure 8). These showed that with 
sufficiently large damping, an absence of resonance was predicted. 


Figure 8 

The predictions of 
Equation (3) for different 
values of the damping 
ratio a. 
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Exercise 2 
(i) Show that the equation of motion for the system in Figure 2 is given by Equation (1). 


(ii) By comparison with the analysis carried out in Section 2, show that the sinusoidal term 
in the steady-state solution of Equation (2) has the amplitude ratio described by 
Equation (3). 

(iii) State a lower bound for the values of the damping ratio a which correspond to 


non-resonant behaviour in this model system. 


[Solution on page 44 | 


Data taken from the apparatus of Figure 6 are shown in Figure 9 below. It was found 
that a reasonable fit to these data could be obtained with the graph for a = 0.2 in 
Figure 8. The prediction of no resonance for sufficiently large damping was also tested. 
With the dashpot moved considerably to the left, a sweep of input frequencies revealed 
a steady decrease of output amplitude with increasing frequency. This indicated an 
absence of resonance, in qualitative agreement with the prediction. 


AlAy 


0 10 20 30 40 Ww 


Figure 9 Amplitude ratio data obtained from the apparatus of Figure 6. 
Part 3 of the programme dealt with the forced vertical oscillations of the tone-arm 
of a record player. It was shown that the large demonstration apparatus of Figure 10 


exhibited damped vibrations when unforced and resonant vibrations when forced. Like 
the previous apparatus, this one involved the pivoting motion of an arm. 


hinged arm 


springy connection 


motor-driven 
’ forcing point 


Figure 10 Demonstration apparatus for Part 3 of the programme. 
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The motion of the tone-arm was modelled by the system shown in Figure 3, in which 
the connection between the stylus and tone-arm is assumed to involve a perfect dashpot 
(linear damping) as well as a perfect spring. In the following exercises you are asked to 
find the equation of motion of this system, and the amplitude of the steady-state 
output when it is subjected to a sinusoidal input vibration. Note that interest centres 
on the motion of the tone-arm relative to the record surface, as described by the 
variable z = x — y. 


Exercise 3 
(i) Show that the equation of motion of the model system shown in Figure 3 is 


mi+re+kxz=klo -—mg+ryt+ ky. 


(ii) Hence show that the displacement z = x — y of the particle relative to the forcing point 
satisfies the differential equation 


mz+rz+kz=klo — mg — my. 


Exercise 4 


If the system shown in Figure 3 is subjected to a sinusoidal input vibration y = Apo cos wt, 
show that the amplitude ratio for steady-state vibrations of the particle relative to the forcing 
point is given by 


Ae (w/o)? : 
Ao \/(1 = (w/o)??? + 402 (w/o)? 


where wo = \/k/m is the undamped angular frequency and a = r /(2V mk) is the damping 
ratio. 


[Solutions on page 44 


In the previous exercise you derived an expression which predicts the amplitude ratio 

for the tone-arm. Some corresponding graphs are shown in F igure 11. It was stated 

that the vertical oscillations of a real tone-arm can be reasonably accounted for by this 

model with a value of about 0.2 for a. This is true. at low frequencies, but there are a 

couple of remarks which need to be made. First, the cartridge detects two sorts of 

motion: relative vertical motion between the stylus and tone-arm, and relative Figure 11 

horizontal motion (the existence of these two signals is what makes stereo reproduction The predictions from the 
; as ; result of Exercise 4. 

possible). Secondly, the arm is not perfectly rigid and can therefore flex in response to 

the motion of the stylus (this is important at frequencies higher than those we consider). 

These two factors mean that the real system has more ‘degrees of freedom’ than have 

been incorporated by modelling the arm as a particle moving in one dimension. 


O 1 2k 4 w]e 


Despite these limitations, the model has certain satisfactory features. It can be shown 
that it predicts resonance for a < 1/\/2 (using the method of Exercise 10 in Section 2, 
page 24), and it suggests that we should increase the damping if we want to reduce the 
amplitude ratio of the system near the natural frequency. However, this natural 
frequency is typically 60 radians per second (that is, 60/(27) ~ 10 oscillations per 
second, or 10 hertz). The frequencies at which vibrations are forced by the undulations 
in the record grooves, producing the music that we hear, are between ten and a 
thousand times greater in value. This led to the question of what could be causing the 
record-player to misbehave so badly. 


The fault turned out to be a bad warp on the record. All records are warped to some 
extent, causing potentially an unwanted resonant response to forcing at low frequencies. 
This effect can be reduced by increasing the damping. It was shown that a similar 
tone-arm on a second turntable did not exhibit such large forced vibrations. The 
difference between the tone-arms was that the second one had a specially designed 
miniature dashpot connecting the arm to the surface of the record. (To make the 
comparison a fair one, a sliding weight on the second arm had been adjusted to 
compensate for the additional weight of the dashpot.) 
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5 - ind of unrt-exererses 


Exercise 1 


A small heavy object is hung from the ceiling by a spring, and it is found that the spring is 
extended by 19.6 cm in the static position. The object is then attached to the floor by two 
identical dashpots, chosen to give a total effective damping ratio a = 1. 


(i) Devise a model for this system, and obtain its equation of motion, giving numerical 


values for the coefficients. Take the value of g to be 9.8ms ”. 


(ii) One of the two dashpots is now removed. What is the new equation of motion? Predict 
the period of the damped vibrations that can now be observed. If the object is released 
from rest with the spring at its natural length, what is the maximum subsequent spring 
extension predicted by the model? | 


Exercise 2 


Suppose that, in addition to the forces acting on the object in Exercise 1(ii), there is an 
externally-produced magnetic force, acting downwards with magnitude 


Po + Pcoswt (Fa > 0), 
Show that the resulting forced vibrations have amplitude 
ae ee (ON 
GQ=p)? +p? 
where p = w/wo and d is the distance by which a constant downward force of magnitude P 
extends the spring. 


Exercise 3 


Figure 1 shows a model of a seismometer, which consists of a particle of mass m suspended 
inside a rigid box by a perfect spring and attached also to a perfect dashpot. The box is given 
a vertical sinusoidal motion of angular frequency w and amplitude Ao, and the displacement z 
of the particle relative to the box is recorded. This displacement is zero when the particle is in 
equilibrium. 


(i) Obtain the equation of motion of the particle relative to the box. 
(ii) Derive the steady-state solution of this equation of motion. 


(iii) Show that, if the damping ratio a and the ratio wo/w of undamped to forcing angular 
frequencies are both small, then z ~ —y. (This result demonstrates that the device can 
be used to record the effect of an earthquake, provided that the damping is light and the 
undamped angular frequency is small compared with the frequencies of ground 
movement. ) 


Exercise 4 


This exercise again refers to the model of Figure 1, and to the results of Exercise 3. VALE ELE bibs etbtitiit 


(i) Show that, if the damping ratio a and the ratio w/wo of forcing to undamped angular y = Ap cos wt 


frequencies are both small, then the relative displacement z of the particle is Figure 1 
approximately proportional to the acceleration of the box, with the constant of 

proportionality being independent of the forcing frequency. (This result demonstrates 

that the device can be used as an accelerometer, provided that the damping is light and 

the undamped angular frequency is large compared with the frequencies of accelerated 

motion to be measured. ) 


(ii) Ifa@=w/wo = 0.1, calculate the ratio of the approximate amplitude of the box 
acceleration (as indicated by the instrument) to the true amplitude. 


[Solutions on page 45 | 


MST204 8 Solutions 


39 


Appendix: Solutions to the exercises 


Solutions to the exercises in Section 1 


Figure 1 


In the equilibrium position, the spring will be extended and 
so the spring force will be upwards. The magnitude of this 
force is T’ = k(l, — Io). The only other force on the particle is 
that due to gravity. As the particle is in equilibrium we have 


T=mq, 
which leads to 
b= lo + mg/k. 


2. The figure required is shown below. The force diagram 
assumes, as stated in the question, that the spring is 
extended and the particle is moving upwards. 


m 


by 


mg 
Figure 2 


As in Exercise 1 we have T = mg when the particle is at its 
equilibrium position. If the particle is below this position 
then, since x < 0, the additional tension in the spring is 
k(—x) = —kz. Alternatively, if the particle is above its 
equilibrium position (but with the spring still extended) 
then, since x > 0, the decrease in tension is kx. In either 
case, the total tension in the spring is given by 


T = mg — kz. 
The particle is moving upwards, so that z > 0 and 
RS Te = se. 
Newton’s second law then gives 
mz =T—R—-—mg 
= —kzx—-rez 
or me +r24+ kz = 0, 
which is Equation (2) of Section 1 once more. 


3. The particle is at a distance x below the equilibrium 

position, which is at a height h above the floor. So the height 

of the particle above the floor is X = h — x. Hence we have 
(= h— XX. ge -X, and f= =X. 

Substituting these expressions into the equation 
mi+rz+kx =0, 

we obtain 


m(—X)+r(-X)+k(h — X) = 
or wmX +X 2X = kh: 


4. ‘The angular frequency 2 of damped oscillations is 


V4mk — r2 ‘ ne 
Q = ——— = 4/u2 - —., 
2m 4m? 
where wo = \/k/m is the undamped angular frequency. 


Hence 2? < w, so that r = 27/Q is greater than the 
undamped period 27/wo. 


5. (i) The damped angular frequency is 
« . Vamk =? 
= 2m 
Here-m = 0.2, k = 2.28 and r = 0.062, so that 


\/4 x 0.2 x 2.28 — (0.062)? 


0.4 
i 
~Y = ~ 3.37rads_t. 


(ii) The undamped angular frequency is 


ck | 2.28 s 
Wo = reese 99 ~ 238rads “g 


6. Equation (7) of Section 1 is 

a(t) = Ae~™ cos(Nt + ¢). 
The derivative of this is 

a(t) = —pAe™ cos(Qt + ¢) — NAe~” sin(Nt + ¢). 
Using the given initial conditions, we have 

xz(0) = Acos¢ = —0.05, (1) 

«(0) = —pAcos¢ — QAsing = 0. 9) 
Now {2 ~ 3.37 from Exercise 5(i) and, using the first of 
Equations (8) in Section 1, 

- 0.062 

i ae = en = Ulu, 

Use of Equation (1) to substitute for Acos ¢ in Equation (2) 


therefore leads to 
Aidice 0.05p oe 0.05 x 0.155 


| ae B.of 
Equations (1) and (3) can be solved by the method of Unit 7 
Subsection 2.3, yielding 

A = \/(—0.05)2 + (0.0023)? ~ 0.050, 

@ = arccos(—0.05/A) ~ 3.1. 


~ 0.0023. (3) 


7. The motion will be oscillatory if r? < 4mk. Here we have 


r=9x10°, m=800, k=5~x 10%, 
so that 
r? = 8.1 x 10’, 


4mk = 4 x 800 x 5 x 10* = 1.6 x 10°. 


Hence the inequality is satisfied, and the motion is 
oscillatory. 
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The ratio of displacements one period apart is 


where 
r 27 V4mk — r2 
p=—, T=— and Q= ——. 
2m Q 2m 
Therefore 
en a 2a x 9 x 10° 
q V4mk — r2 /1.6 x 108 — 8.1 x 10° 
4 
eo hae 
8.89 x 10° 


From this it follows that 
eee ise ee ee 


Hence one complete cycle after being released the 
displacement of the car body below its rest position will be 


0.1 < 1.73 « 10° =1.73-x 10° “m, 
or about 0.2 mm. 


8. For the differential equation 
mz+rz+kx =0, 

the auxiliary equation is 
mr +rA+k=0 

with roots 


oe —r+t+/r2 —4mk 
aay a 


In each of the following three cases we apply Procedure 1.1 
of Unit 6. 
(i) 1? < 4mk (complex roots for \) 
The roots are \ = —p +72, where 
i V4mk — r? 
 seciite eis Wee 
2m 2m 

so the general solution of the differential equation is 

a(t) = Ae” cos(Qt + ¢). 
Here A and ¢ are arbitrary constants, though all solutions of 


this form may be obtained with A non-negative and ¢ 
between —7z and z. This case is called weak damping. 


and t= 


y) 


(ii) r? = 4mk (equal roots for 4) 
The auxiliary equation has just one solution, \ = —p, where 
p=r/(2m). So the differential equation has general solution 
ails Be" +O. 
where B and C are arbitrary constants. This case is called 
critical damping. 
(iii) r? > 4mk (distinct real roots for 2) 
The roots are \ = —p; and A = —p2, where 
r— Vr? —4mk r+w~J/r?2—4mk 
2m 2m 
Note that both p; and p2 are positive. The general solution is 
Pie. ee, 
where B and C are arbitrary constants. This case is called 
strong damping. 


pi = and p2= 


9. (i) The most efficient way of finding the stiffness k is 
to use energy conservation (Unit 7 Section 3), since there is 
no damping in operation between the firing of the gun and 
maximum spring compression. Equating the initial kinetic 
energy, with m = 10° and uv = 30, to the potential energy for 
a compression of 2 metres, we obtain 


7 eID x SU = ee 

ar k= 2965 Na; 

(ii) For critical damping we require r? = 4mk, that is, 
r = 2V'mk = 21/108 x 2.25 x 105 = 3 x 10°. 


Thus a damping constant of 3 x 10* Nm7's is required. 
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Solutions to the exercises in Section 2 


1. We look for a particular (steady-state) solution of the 
form 


x = Bcoswt+ Csinuwt, 
whose first two derivatives are 

z= —Bwsinwt + Cw cos wt, 

7% = —Bw’ coswt — Cw” sinwt. 
Substitution of these expressions into the differential 
equation 

mz+rxe+kx = Pcoswt 
produces 

m(—Bw* coswt — Cw” sin wt) 

+ r(Cw cos wt — Bw sin wt) 

+k(Bcoswt + Csinwt) = Pcoswt. 


We can now compare the coefficients of the cosine and sine 
terms in turn on both sides of this equation: the coefficient 
of coswt on the left-hand side must equal P on the right, and 
the coefficient of sinwt on the left-hand side must equal zero 
(since there is no sinwt term on the right). This gives 


—wee +#Cu + ke = FP, (1) 
Cw -~ Tie =f. (2) 
From Equation (2) we have (k — mw*)C = rwB, or 
rwB | 
ae, a 
2 k — mw? (3) 
Substituting for C in Equation (1), we obtain 
2 r?w? 
(k — mw") B+ 7— 5B=P 
which leads to 
P(k — mw?) 


Se 


It follows from Equation (3) that 
Prw 


Having found the steady-state solution in the form 


C= 


x = Bcoswt+ Csinuwt, 
it remains to translate this into the form 
x = Acos(wt + ¢), 
where, from Unit 7 Subsection 2.3, 
A= 7B? 4C-. 
and = -—arccos(B/A) (since C > 0). 


Using the expressions above for B and C, we find 


ee 2 
and = -—arccos ar Sal oe ; 
(k — mw?)? + r2w? 


in agreement with Equations (3) and (4) of Section 1. 


[Notice that this approach takes longer than the phasor 
method. | 


2. (i) The angular frequency is 
w=2xf (where f is frequency in Hz) 
=2nx2.5~15.7rads’. 


(ii) From Equation (3) of Section 2, the amplitude is 
ee ? 
With m = 800, k =5 x 10*, r = 9 x 10°, P = 800 and 
w ~ 15.7, we have 
k — mw? ~ —1.47 x 10° 
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and 
> 800 


~ 4/(=1.47 x 10°)? + (9 x 108 x 15.7)? 


800 800 es 
Y ———  ~ ———_ ~ 3.92 x 10 
417x101 2.04 105 - m, 


so the amplitude of the vibration is about 4mm. 


(iii) From Equation (4) of Section 2, the phase angle of the 


vibration is 
k — mw? 
(k — mw?)? + r2w? 


o@ = — arccos ( 


Using the values above, we obtain 
@ ~ — arccos e. : =| 
2.04 x 10° 
~ — arccos(—0.722) 
~ —2.38rad (or —136°). 
The time lag is therefore 


3. The applied force is cost = Re(e”), with angular 
frequency w = 1. We look for a steady-state solution of the 
form x = Re(Ze"'), where Z is a complex constant (namely, 
the phasor of x). Successive differentiations of this give 


& = Re(iZe") and = Re(—Ze"’). 
Substitution of these expressions into the equation of motion 
produces . 


Re(—2Ze” + 5iZe”* + 14Ze") = Re(e’) 
or Re((12+5i)Ze) = Re(e”’). 
This equation will be satisfied provided that 


(12 + 52)Z = 1, 
for which we require 
1 SL eee oe 
{fi eS 
where 
es 12 
¢@ = — arccos(35) 


~ —0.395rad (or —22.6°). 
Hence the steady-state solution is 
xz = Re (e'%e") = Re (Ae) 
= 4+ cos(t + 4), 
where @ ~ —0.395 rad. 


4. (i) Using Euler’s formula to expand e’“*, we have 
(M —iN)e*™ 
= '(M —iN)(coswt + isinwt) 
= (M coswt + Nsinwt) + i(M sinwt — N coswt), 
so that 
Re((M —iN)e**’) = Mcoswt + N sinwt. 
(ii) The applied force is (using the above) 
2 cos 2t + 9sin 2t = Re((2 — 9i)e?”*). 
The angular frequency of this applied force is w = 2, so we 


look for a steady-state solution of the form 2 = Re(Ze”"*), 
whose first two derivatives are 


& = Re(2iZe*") and # = Re(—4Ze7"). 
Substituting these expressions into the equation of motion, 
we have 
Re(—4Ze™™ + 4iZe” + 5Ze7"*) = Re((2 — 9i)e””) 
or Re((1+4i)Ze?”) = Re((2 — 9i)e?"). 
This equation will be satisfied if 
(1+ 41)Z = (2 — 9), 
for which we require 
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ie Ok Ae At 

eg ae 12 + 4? = : 
Using the result of part (i) once more, the steady-state 
solution is 


x = Re((—2 — i)e*"*) 
= —2cos 2t + sin 2t. 


To find the amplitude A and phase ¢ of this vibration, we 
express the phasor Z in its exponential form Z = Ae’?. This 


gives 
A= \/(—2)2? +1? = V5 ~ 2.24m 
and = -—arccos(—2/V5) ~ —2.68 rad, 


so the amplitude and phase are respectively 2.24m and 
—2.68 rad. 


<2 4, 


5. We seek the phasor Z for the steady-state solution of 
P 
%+2awot + usr = —coswt = Re (—e'“*) 
m m 
Putting ¢ = Re(Ze““") leads to « = Re(iwZe***) and 
&% = Re(—w? Ze"), so the phasor Z is given by 
(—w? + i2awow + wo)Z = P/m 
os 
wg — Ww? + i2awow 
— P(wp — w? — i2aww) 
~— m((w2 — w?)? + 4a2w2w?) ’ 
Hence the amplitude of the steady-state solution is 
Py/(w2 — w?)? + 4a2wew? 
m((w2 — w?)? + 4a2wew?) 
P/m 
a/ (wa — w?)? + 4a2wew? 


and the phase is 


A= = 


@ = Arg(Z) = — arccos ( 


Wo — iar 
(we — w?)? + 4da2wew? 


6. Dividing each coefficient in the differential equation by 
the coefficient of @ (so that the coefficient of z becomes 
unity), we obtain 


E+ 2z + 252 = 2 cos 6t. 
We now compare the coefficients of this equation with those 
of Equation (10) in Section 2, from which it follows that 
3, we = 25, P/m= 2. 
We also deduce from the two right-hand sides that w = 6. 


(i) From above, the undamped angular frequency is wo = 5, 
so that the damping ratio is 
9 9 
= — = — = 0.225. 
ee 
(ii) The amplitude of the steady-state solution is 
P/m 
(we — w?)? + 4arwew? 
2.25 
(25 — 36)2 + (2.25 x 6)? 
2.25 


= ————_—_—. ~ 0.129 m, 
V121 + 182.25 


and the phase is 


ee 
wo — wW 
6 = <areces | 
(we — w?)? + 4arwew? 


= — arccos SNe eee 
( sas + a) 
—2.255rad (or —129.2°). 


The angular frequency of the steady-state solution is 
w= 6rads '. 


2awo — 


I? 
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7. From Solution 6 above, we have 
Scenes k 


—=-, — =u = 25, A ~ 0.129. 
m 4 m 
‘Hence k/P = 100/9 and 


r 9 


8. The relative error will be less than 10% provided that M 
is between 0.9 and 1.1, which from the graph for a = 0.6 
appears to be the case for the approximate range 

0<w/wo < 0.9. Since 0.9W9 = 0.9 x 20 = 18, the 
corresponding frequency range (for w) is between 0 and 
18rads~°. 


9. We use the values m = 800, k = 5 x 10*, r = 9 x 10° 
_and P = 800 given in Exercise 2, together with the value 
w ~ 15.7 obtained in Solution 2(i). The undamped angular 
frequency is 


rk /5 x 104 i 
Wo = ag os 300 ~ 7.91rads es 


and the damping ratio is 
eee 9 x 10° 
2/mk 2/800 x 5 x 104 


[Note that this value of a is just above 1/\/2, so that there is 
no input frequency which makes the car body resonate.| 


ants US PT 


The relevant value of w/wo is approximately ee 28. The 
corresponding value of M from Figure 5(a) in Section 2 is 
approximately 0.25. Since M = Ak/P, we then have 
MP _ 0.25 x 800 
ko eee 
The estimated value of ¢ from Figure 5(b) is 

o@~ —0.75r ~ —2.36 rad. 
The estimated values of A and ¢ here are both close to the 
values calculated in Exercise 2. 


a —4x 107°. 


10. With p = w/wo, Equation (13) of Section 2 becomes 
1 
M-= 


/(1 — p?2)2 + 4a2p2 
Differentiation with respect to p gives 
dM _ —3[2(1 — p?)(-2p) + 80”p] 


dp [(1—p?)? + 4a2p?]8/2 
so that dM/dp = 0 (for points on the graph with horizontal 
slope) when 


— 4p(1 — p*) + 8a*p = 0 
or p(p* —1+ 2a”) =0. 
The roots of this cubic equation are 


p= 0 -and- p= 44/1 — 2a2. 


(i) We know that p = w/wo is non-negative. For a < 1/V2, 
the non-negative roots are p = 0 and p= v1 — 2a?, and the 
graphs in Figure 5(a) of Section 2 indicate that the second of 
these gives a maximum value for M. (This can be confirmed 
by examining the sign of dM/dp to the left and to the right 


of p = V1 — 2a?.) 

(ii) For a = 1/V2, the root p = 0 is repeated, and for 
a> 1 V2, the expression 1 — 2a? gives complex roots. 
Figure 5(a) shows that in this case M is a decreasing 


function of p for p > 0, so that p = 0 gives the maximum 
value for M. (It can also be shown from above that 


dM /dp < 0 for p> 0 and a > 1/v2.) 


MS7T204 8 Solutions 


Solutions to the exercises in Section 3° 


Figure 1 


In deriving the equation of motion we assume that the spring 
is extended and that the dashpot’s length is increasing. The 
corresponding force directions are as shown above. 


The extension of the spring is x — /p, and the rate of increase 
of the dashpot’s length is «. The spring and dashpot forces 
therefore have respective magnitudes 


f ere k(x ae lo) 
The only other force acting on the particle is the 


gravitational force, of magnitude mg vertically downwards. 
Newton’s second law therefore gives 


Sak. = og 


mz=mg—T—R 
= mg — k(x -—I9) — rz 
or mz+r¢z+krx=mg-+ klo. 
This equation of motion is identical to Equation (3) in 


Section 1, as expected (since the model considered is the 
same in the two cases). 


Figure 2 


We assume, as usual, that the spring is in tension and that 
the dashpot is lengthening. The corresponding force diagram 
is shown above. Note that Pcoswt is the x-component of the 
external force, so that this force has the direction shown only 
when P cos wt is positive (the same applies to the diagram 
given in the question). 
The spring force has magnitude T = k(x — Io) and the 
dashpot force has magnitude R = rz, so that the required 
equation of motion is 

mz = Pcoswt — rx — k(x — Io) 
or mz+re+kzr = klo + Pcosut. 


3. We seek the phasor Z = Ae’®, where x = Re(Ze"’) is 
the steady-state solution of the differential equation 


mz+rz+kzx = Ao(kcoswt — rwsinwt). 
Substituting for x and its derivatives, this equation can be 
written as 

Re[(—mw” + irw + k)Ze"] = Re[Ao(k + irw)e**], 
where the result of Exercise 4(i) of Section 2 has been used 
on the right-hand side. This equation will be satisfied if 


(—mw* + irw +k)Z = Ao(k + irw) 
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or 
Ao(k + irw) 
Fe 
(k — mw?) + irw 
_ Ao(k + irw)[(k — mw?) — irw] 
aS (k — mw?)? + r2w2 
Ao ([k(k — mw”) + r2w?] — imrw®) 
(k — mw?)? + r2w? , 
This is of the form Z = u — iv, where 
_ Ao[k(k — mw?) + r?w?] 
(k= mw?)? + r2w? 


and 
Aomrw® 
(k — mw?)? + r2w2’ 


so that the amplitude is 
A= |2Z| = Vu? + v? 
Ao v/ [K(k — mw?) + r2w?]2 + (mrw3)? 


(k — mw?)? + r2w2 


VU 


and the phase is 
@ = — arccos(u/|Z]) 


= — arccos ( 


4. (i) When the particle is in equilibrium, the dashpot 

length is constant and so the dashpot exerts no force. The 
force diagram is therefore as shown below, where T is the 

magnitude of the spring force. 


7 


k(k — mw?) + r?w? 
[A(k — mw?) + r2w?]2 4 (mrw?)? 


£ 


Figure 3 


In equilibrium these two forces must balance, so that 
I’ = mg. Now if the spring has length /; when the particle is 
in equilibrium then its extension is J, — lo, giving 
T = k(li — lo). Hence 
k(1y = lo} = mg, 
from which the equilibrium length of the spring is 
li =Ilo + mg/k. 
(This is the same argument as that required to answer 
Exercise 1 of Section 1.) 


(ii) : 


Figure 4 (a) i, 
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Figure 4(a) shows the system when displaced from its rest 
position. The displacements x of the particle and y of the 
forcing point are as described in the question. F igure 4(b) 
shows the forces which act on the particle when the spring is 
extended (as it must always be if a rope is being modelled) 
and the particle is moving upwards (that is, the dashpot is 
contracting). The magnitudes of the spring and dashpot 
forces are denoted by T and R respectively. 
The extension of the spring is 

lo+mg/k+y—2x-—lo, 
so that 

T= k(mg/k+y-— 2). 
Since © > 0, we also have R = r&. Application of Newton’s 
second law gives 

mz =T— R—mg 

= hlimg/k +y 2) — re — ma, 

OF me re 4 ka = ky: 
This is the required equation of motion. 
(iii) Putting y = Ao coswt in the result of part (ii) gives the 
equation of motion 

mz+rz+kx = kAo coswt. 
This is the same as Equation (1) of Section 2, except for the 
presence of kAo rather than P on the right-hand side. The 
graphs in Figure 5 of Section 2 can therefore be applied to 
this situation where, with k Apo in place of P, we have 
M = A/Ao. Note that M here is the ratio of the 
displacement amplitudes of the particle (output) and of the 
upper end of the spring (input). (This quantity is called the 
amplitude ratio of the system.) 


From the values given in the question, we obtain 


wo = /k/m = V50 ~ 7.07 rads ', 
w=2x 27 ~12.6rads’, 


and hence w/wo ~ 1.78. Corresponding to this value of w /wo, 
the graph for a = 0.2 in Figure 5(a) of Section 2 gives 
M ~ 0.45, from which 

A =~ 0.45Ao = 0.45 x 0.05 ~ 0.023. 
The amplitude of the particle’s steady-state oscillations is 
therefore about 0.023 m, or 2.3cm. 


[We may check that the spring representing the rope is in 
tension throughout the motion by showing that T (as defined 
above) is always positive. Since |z| < A and |y| < Ao, we 
have | 

l= mg + k(y — x) > mg — k(Ao + A). 
The value of the right-hand expression is approximately 

4 x 9.8 — 200(0.05 + 0.023) ~ 24.6N, 


so the spring is always in tension.] 
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Solutions to the exercises in Section 4 


wai 


ee: ed ar a 
oe ee 


ee 1 2 2 4 
Figure 1 


The model predicts that a graph of 

log.{(Ao + Ao)/(An + A;,)] against n should be a straight 
line through the origin (with slope pr). The given data have 
been plotted on such a graph above. There is no straight line 
through the origin which provides a convincing fit to these 
data points, so this test indicates that linear damping is not a 
satisfactory modelling assumption for a water-filled dashpot. 


2. (i) Assume that the spring is extended and that the 
dashpot is lengthening. The spring force is then of magnitude 
T = k(x — y — Io) and the dashpot force is of magnitude 
R=rz, both being directed downwards. The only other 
force acting is that due to gravity, which is of magnitude mg 
directed downwards. Newton’s second law therefore gives 

me = —T Kh =mg 

= —k(x4 —y—Io) —rz — mag, 

or mz+rz+kz=klo —mg+ky, 
which is Equation (1) of the text. 
(ii) Equation (2) of the text is 

mz+rz+kxz = k(lo + yo —mg/k) + kAo cos wt. 
The steady-state solution is 

x =I9 + yo — mg/k + Acos(wt + 9), 
where the sinusoidal term is the steady-state solution of the 
equation 

mz+rz+ kx = kAocoswt. 
We now proceed much as in Solution 4(iii) of Section 3. The 
last equation is the same as Equation (1) of Section 2, with P 


replaced by kAo. The analysis of Subsection 2.3 can therefore 
be applied, leading to Equation (13) of Section 2, which was 


1 
(1 — (w/wo)?)? + 40?(w/wo)? 
Here M is equal to A/Ao, since P has been replaced by k Ao. 


This is the expression for the amplitude ratio given in 
Equation (3) of Section 4. 


(iii) From the result of Exercise 10 of Section 2, there will be 
no resonance provided that the damping ratio a is greater 
than or equal to 1/,/2. 
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eas ; Oe 


Figure 2 


We assume that the spring is extended and that the dashpot 
is lengthening. Then the spring and dashpot forces are as 
indicated in the diagram. The lengths of the spring and 
dashpot are both x — y. Hence the extension of the spring is 
x — y — lo, and the rate of increase of the dashpot’s length is 
x —y. The magnitudes of the spring and dashpot forces are 
therefore 


T=k(a@—y-—Io) and R=r(e-y) 
respectively. So the equation of motion of the particle is 
mz = —-T—-R—-mg 
= —k(x —y— lo) —r(z— y) — mg 
or mz+ra+kr=klo —mgt+ryt+ky. 
(ii) Since z = x — y, we have 
¢=2—9 and. 2=2— 4%. 
The equation of motion becomes 
mzZ+y)+r(ié+y)+k(et+y) =klo-—mgtryt+ky 
or mz+rz+kz=klo —mg —my. 


4. The equation of motion of the tone-arm in terms of 
z=a2-—y is given by the last equation in Solution 3 above. 
With y = Ao coswt we have ij = —w” Ao coswt, giving 
mi+tritkz = klo — mg + mw” Ao cos wt. 
The steady-state solution of this equation is 
z =I —mg/k + Acos(wt + ), 
where the sinusoidal term is the steady-state solution of 
mi +ri+tkz = mw’ Apo cos wt. (1) 
This is, once again, the same as Equation (1) of Section 2, 


but with mw” Apo in place of P. As before we may therefore 
conclude that Equation (13) of Section 2 holds, which is 


Dl is src Poets at 
(1 — (w/wo)?)? + 40? (w/wo)? 
Here M = Ak/P (from Section 2) is replaced by 
Me sh 

mw? Ao w? Ao : 
since Wo = ,/k/m. Hence we have the amplitude ratio 

tes (w/wo)? 

Ao (= (w/o)? )? + 402 (w/o)? 
as required. 


re 


As an alternative, we show below how the phasor method 
can be used to derive this expression directly, without 
reference to Section 2. Starting from Equation (1) above, we 
divide through by m and then re-express the result in terms 


of the undamped angular frequency wo = ,/k/m and the 
damping ratio a = r/(2Vmk). This gives 

Z+ 2awopz + woz r= we As cos wt. 
We seek a solution to this equation of the form 
z = Re(Ze’”'), where Z is a complex constant. Since 
z = Re(iwZe™") and % = Re(—w” Ze*”"), substitution into 
the differential equation produces 


Re[(—w* + 2iawow + wo)Ze*] = Re(w* Ave”). 
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This is satisfied provided that 

(—w? + 2iawow + we)Z = w? Ao, 
for which we require 

w? Ao 
we — w? + 2iawow 

= (w/wo)* Ao 
1 — (w/wo)? + 2iaw/wo 

The amplitude of these steady-state vibrations is 
(w/wo)* Ao 


V(1 = (w/wo)?)? + 402(w/wo)? 


giving the required amplitude ratio as before. 


 — 


A=|2Z|= 


Solutions to the exercises in Section 5 


s ee it | We model the object by a particle of mass m, and 
assume that the spring and dashpots are perfect, with the 
spring having stiffness k and each dashpot having dashpot 
constant r. The particle displacement x is measured 
downwards from the equilibrium position. This model 
system is shown on the left below. 


7 


‘\ 
R R 
mg 
oo 
Figure 1 | 


Assuming that the spring is extended and the dashpots are 
contracting, the forces acting on the particle are as shown on 
the right of the diagram above. The magnitudes of the 
spring and dashpot forces are respectively 


T=mg+ker 
(since T = mg when the particle is in equilibrium) and 
= Ts = ra. 
Newton’s second law then gives 
mz=mg—T—-R-R 
= mg —(mg+ ka) — 2rz¢ 
or mz+2rz+kzr=0. 
This can be written as 
% + 2awot + wor = 0, 
where wo = \/k/m is the undamped angular frequency and 


a = r/Vmk is the total effective damping ratio (that is, the 
ratio obtained by regarding all the damping as being 
produced by a single dashpot with dashpot constant 27). In 
the equilibrium position, a force of magnitude mg produces a 
spring extension of 19.6cm = 0.196 m, so that 


mg = 0.196k 
ae eae ae 
CS ae 


Also a = 1 is given, so the equation of motion is 
£ + 2V 502 + 50z = 0. 
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(ii) Removing one dashpot has the effect of halving the 
damping. The new equation of motion is therefore 


+ /50¢ + 50zr = 0, 


where the damping ratio is now a = r/(2Vmk) = 0.5. From 
Equations (7) and (8) of Section 1, the general solution of 
this equation is 

a = Ae~™ cos(Nt + ¢), (1) 
where 


= 2 am 7" 


2m 2m 


and A, ¢ are arbitrary constants. From w% = 50 and a = 0.5, 
we have 


k= 50m ~and r = Vmk = V50m, 
giving 

p=V12.5 and 2 = 4200-50 = V3755. (2) 
The period of damped vibrations is 


= wi ona oT 03's. 


hae 
If the particle starts from rest with the spring at its natural 
length then, since the origin is at the particle’s equilibrium 
position, the initial conditions are 
z(0) = —0.196 and 2z(0)=0. 
Putting ¢ = 0 into Equation (1) and applying the first of 


these conditions produces 


Acos @ = —0.196. 
The maximum subsequent spring extension occurs after half 


a period at time t = $7 = 7/Q, when x(t) = 0 once more. 


The corresponding value of x is (using Equations (2)) 
a(r/Q) = Ae~°*!™ cos(x + ¢) 
= — Ae °*/® cog o) 
= 0.196e~7/¥3 ~ 0.032. 
The spring extension is then approximately 
0.196 + 0.032 = 0.228 m, 
or approximately 22.8 cm. 


2. The equation of motion for Exercise 1(ii) was 
mz+re+kx =0, 
and with the additional downward magnetic force, this 
becomes 
mz+rz+kxz = Po + Pcosut. 
The steady-state solution is 
x = Po/k + Acos(wt + ¢), 


where the sinusoidal term is the steady-state solution of the 
equation 


mz+ra2+kxz = Pcoswt, 
that is, 


ts 
# + 2awot + wox = — cosut, (3) 
m 


where wo = \/k/m and a = r/(2Vmk). Equation (3) is 


identical to Equation (10) of Section 2, and the 
corresponding amplitude of forced vibrations is given by 
Equation (11) of Section 2 as 


P/m 
(w2 — w2)? + 4a2wew? 
Now a = 0.5 from Exercise 1(ii), so the amplitude is 
P/(mw9) 
(1 — (w/wo)?)? + (w/wo)? 
eles Gee 
where p = w/wo and d = P/(muw) = P/k is the distance by 


which a downward force of constant magnitude P extends 
the spring. 
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3. (i) Let x be the upward displacement of the particle, 
measured from the fixed point about which the sinusoidal 
motion y = Ao coswt of the base of the box takes place. 
Then we have x = y + z +d, where d is the length of the 
dashpot in the equilibrium position. Assuming that the 
spring is extended and the dashpot lengthening, the forces 
acting on the particle are as shown in the diagram below. 


Figure 2 


The spring and dashpot forces are respectively 

i = ig — kz 
(since J’ = mg when the particle is in equilibrium) and 

i = Fe. 
Newton’s second law therefore gives 

mz =T—mg—R 

= (mg — kz) — mg — r. 

On putting c = y+ z2+d we have ¢ = y + 2, so the equation 
of motion relative to the box is 

mz+trz+kz=—my 
or mi+ri+kz = mw’ Ao cosut. 
(ii) The last equation is a copy of Equation (1) of Section 2, 
with mw” Ao in place of P, so by Equations (11) and (12) of 
that section we have the steady-state solution 

z= Acos(wt + 9), 
where 
igs 


(we — w?)? + 4a2wew? 


2 2 
and ¢=-—arccos | —— 
(we — w?)? + 4a2wew? 


A= 
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(iii) If wo/w and a are both small, then 
A ey, ' 
((wo/w)? — 1)? + 4a?(wo/w)? % 


and 


2 @® Beem 


@ = — arccos (ete 
((wo/w)? — 1)? + 4a? (wo /w)? 
~ —arccos(—1) = —7. 
Hence we have 
z = Acos(wt + ¢) 
~ Ao cos(wt — 7) = —Ap coswt 


or 2 = yy. 


4. (i) Ifw/wo and a are both small, then 
(w/wo)* Ao is w* Ao 


2 


A= es 
(1 — (w/wo)?)? + 40?(w/wo)? = a 


and 


@ = — arccos (ott | 
(1 — (w/wo)?)? + 40?(w/wo)? 
~ —arccos(1) = 0. 
Hence we have 
z= Acos(wt + ¢) 
~ (w/wo)* Ao cos wt 
or z 2 —Yy/w9. 
Thus z is proportional to y, and the constant of 
proportionality —1/wé is independent of w. 
(ii) The approximate value of acceleration as indicated by 
the instrument is 
Wee = we cos(wt + ¢), 
whose amplitude is wo A. The actual value of the acceleration 
is 
yy = —w’ Ao cos wt, 


with amplitude w? Ao. Using the expression above for A and 
putting a = w/wo = 0.1, the ratio of these amplitudes is 


ar = 1 

wAo — \/(1 — (w/wo)?)? + 402(w/wo)? 
evs 1 
~ ./— (0.12)? +4 x 0.14 
one ii fe 


oe 
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